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PREFACE 
This  volume  is  the  third  and  last  of  a  series  reports  of  a  project  funded 
by  the  Environmental  Protection  Agency.   An  overview  of  the  project  objectives 
and  accomplishments  is  given  in  Chapter  One  in  Volume  One.   In  the  present 
volume,  the  first  chapter  was  written  by  L.  0.  Myrup,  the  second  by  L.  0. 
Myrup  and  J.  Lague  and  the  third  by  L.  0.  Myrup,  K.  Lidster,  D.  C.  Hudson  and 
W.  B.  Goddard.   Publication  of  this  report  has  been  delayed  for  over  a  year 
in  order  that  verification  of  the  line  source  model  discussed  in  Chapter  Three 
could  be  included.   Unfortunately,  this  was  not  possible  due  to  the  complexity 
of  the  appropriate  experiments.   Therefore,  Chapter  Three  consists  only  of  a 
description  of  the  development  of  the  model.   We  plan  to  include  verification 
information  when  the  model  is  published  in  the  literature. 
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CHAPTER  ONE 

"THE  IMPACT  OF  HIGHWAYS 
ON  AIR  QUALITY" 

L .  0 .  Myrup 
August  1972 


Background 

The  twentieth  century  has  been  a  time  of  constantly  increasing  energy  use 
in  the  United  States  and  throughout  the  world.   In  the  United  States  since 
1925  the  rate  of  increase  has  been  about  3%  per  year  (Darmstadter ,  1971).   At 
the  same  time  there  has  been  a  major  shift  from  coal  to  petroleum  products  as 
the  major  source  of  energy  in  this  country  and  throughout  the  world.   Pollution 
is  inherent  in  the  production  and  utilization  of  energy.   Ultimately,  we  must 
bow  to  the  Second  Law  of  Thermodynamics  and  recognize  that  waste  heat  is  the 
cleanest  kind  of  pollution  we  can  possibly  hope  to  achieve.   Of  course,  the 
relative  efficiency  of  various  kinds  of  energy  production  and  use  varies  and 
is  directly  related  to  resultant  pollution.   In  general  terms,  the  most  inef- 
ficient uses  which  we  make  of  energy  resources  are  (1)  the  conversion  process 
to  electric  power  and  (2)  transportation  by  means  of  the  internal  combustion 
engine.   Both  in  terms  of  waste  heat  and  pollution,  power  generation  and 
transportation  lead  all  other  categories  (Darmstadter,  1971;  List,  1971). 

The  subject  of  this  paper  is  pollution  associated  with  transportation,  in 
particular  that  emitted  by  motor  vehicles  traveling  on  highways.   In  the  course 
of  the  discussion  we  shall  also  consider  the  relative  importance  of  pollution 
emitted  by  residential  and  arterial  streets  in  addition  to  highways.   The 
processes  which  lead  to  a  specified  air  quality  level  in  the  air  which  we  all 
breathe  is  complex.   Figure  1-1  is  a  rough  sketch  of  the  kind  of  perspective 
which  is  useful  in  considering  this  problem.   The  diagram  suggests,  and  we 
are  all  aware,  that  the  meteorological  influences  are  highly  important  in  de- 
termining air  quality.   I  have  chosen  to  emphasize  in  this  paper  the  meteoro- 
logical aspect  of  the  problem,  particularly  the  current  state  of  our  ability 
to  model  and  predict  the  pollution  emission   to   air-quality  part  of  the 
process.   In  addition,  I  will  discuss  the  influence  of  land-use  patterns  on 
the  microclimate  and  pollution  emission.   I  will  not  attempt  to  deal  with  the 
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Figure  1-1.  The  air  pollution  systan. 


human  health  problem  associated  with  air  pollution. 

Much  of  this  paper  is  meant  to  have  implications  for  and  to  be  useful  to 
the  policy  decision  and  planning  processes  in  society.   One  of  the  benefits 
of  making  a  diagram  like  Figure  1-1  is  that  it  forces  you  to  consider  the 
problem  from  a  broad  perspective.   From  this  point  of  view,  it  becomes  clear 
that,  say,  the  inventory  process  or  land-use  patterns  are  also  critical  com- 
ponents in  determining  the  net  result  of  the  entire  process  which  I  take  as 
human  health  and  comfort.   The  policy-decision-making  and  planning  processes 
are  particularly  important,  for  it  is  there  that  we  have  some  hope  that  ac- 
curate information  and  our  best  estimates  of  the  consequences  of  changes  in 
specific  portions  of  the  system  may  have  significant  impact  for  the  common 
good. 

In  this  context,  I  would  like  to  point  out  that  the  transportation  in- 
dustry has  a  unique  advantage  and  role  to  play.   In  contrast  to  many  portions 
of  society,  the  advance  planning  component  of  transportation  is  strong  and 
accustomed  to  making  us  of  the  best  available  engineering  information  in  its 
operations.   Consequently,  the  additional  considerations  which  are  necessary 
to  protect  and  manage  air  quality  come  naturally,  in  my  opinion,  to  highway 
engineers  and  administrators.   My  optimistic  view  is  then  that  the  transporta- 
tion industry  has  an  excellent  opportunity  to  lead  the  way  in  rational  plan- 
ning to  simultaneously  satisfy  both  the  innemse  demands  which  society  makes 
for  goods  and  services  and  also  to  protect  environmental  quality. 
Discussion  of  the  Problem 

Since  we  are  dealing  with  pollutants  emitted  from  motor  vehicles,  we  may 
technically  define  the  problem  from  a  meteorological  point  of  view  as  trans- 
port and  diffusion  of  material  emitted  from  a  line  source  near  the  surface  of 
the  earth.   The  problem  is  not  really  this  simple,  as  we  shall  see,  but  the 
geometry  strongly  suggests  this  approach.   One  of  the  problems  is  that  of 


TABLE  1-1. 
Highway  height  %Improvement  per  meter  height 

Z(m)  "  10°  9C   mo/   i    / 

— —  —  =  100/z  In  z/z 

L      <3z  o 


1  43.5 

2  16.7 

3  9.8 

4  6.8 

5  5.1 

6  4.1 

7  3.4 

8  2.9 

9  2.5 
10  2.2 


scale.   Like  most  complex  systems,  the  factors  which  are  important  to  small- 
scale  processes  in  the  atmosphere  are  very  different  from  those  which  de- 
termine the  larger  scales.   The  phenomena  which  effect  the  transport  and  dis- 
persion of  highway  pollution  include  the  micro-  and  mesoscale,  in  the  the 
meteorological  terminology.   In  this  section,  we  will  develop  this  distinction 
and  its  significance  to  highway  air  quality. 

Eschenroeder  (1971)  has  shown  the  existence  of  a  zone  above  highways  in 
which  the  air  is  well  mixed  by  the  energy  of  the  moving  vehicles.   This  zone 
has  the  dimensions  of  twice  the  height  of  the  average  vehicle  times  the  high- 
way width  and  is  called  the  mechanically  mixed  cell.   This  mixing  zone  or  cell 
can  be  thought  of  as  the  cumulative  effect  of  all  the  turbulent  wakes  which 
are  formed  behind  each  object  moving  along  the  highway.   Within  the  mechanically 
mixed  cell,  pollution  is  dispersed  uniformly,  to  a  first  approximation,  and  this 
state  is  the  effective  initial  distribution  of  air  pollution  which  is  acted 
upon  by  micro-  and  mesometeorological  processes  to  affect  air  quality  down- 
wind of  the  highway. 

It  should  be  noted  that  in  this  simple  model  the  quality  of  the  air  which 
drivers  and  passengers  breathe  directly  above  the  highway  is  completely  de- 
termined by  the  properties  of  the  vehicles  and  their  motors  and  is  independ- 
ent of  meteorological  processes.   At  sufficiently  high  winds,  this  certainly 
cannot  be  true.   The  extent  to  which  the  mechanically  mixed  cell  is  affected 
by  external  meteorological  influences  has  not  been  studied  in  sufficient  de- 
tail to  allow  any  general  statements  at  this  time.   This  is  an  important  point 
because  it  determines  the  pollution  levels  which  highway  travelers  are  sub- 
jected to  and  should  be  the  object  of  an  intensive  experimental  effort. 

One  operational  conclusion  which  can  be  drawn  from  the  mixing  cell 
concept  is  that  the  effective  height  of  emission  of  highway  pollutants  is  the 
average  height  of  the  vehicles.   For  the  remainder  of  this  section,  we  shall 
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imagine  that  the  actual  highway  emission  pattern  may  be  replaced  by  a  line 
source  concentrated  at  that  height.   This  material,  whether  gaseous  or  par- 
ticulate, is  now  subjected  to  the  microscale  wind  transport  and  diffusion 
processes.   Under  most  conditions,  the  wind  near  the  ground  varies  logarith- 
mically with  height  according  to 

U* 

u  =  r  mf-  ,d) 

o 

where  U  is  the  wind  speed  at  height  z ,  k  is  the  von  Karman  constant,  U.  is  a 

micrometeorological  parameter  called  the  friction  velocity,  and  z   is  the 

o 

aerodynamic  roughness  and  is  directly  related  to  the  geometrical  roughness  of 
the  surface.   These  last  two  parameters  are  important  to  the  processes  by 
which  atmospheric  turbulence  is  generated  near  the  ground.   U. ,  which  is  de- 
fined by 

where  x  is  drag  force  which  the  wind  exerts  on  the  surface  and  p  is  air  den- 
sity, can  be  thought  of  as  representative  of  the  turbulent  velocity  fluctu- 
ations which  arise  from  the  roughness  of  the  surface  and  which  do  most  of  the 
work  in  diffusing  material  near  the  surface.   The  significance  of  the  logar- 
ithmic wind  law  for  highway  air  quality  lies  in  the  fact  that  in  all  theo- 
retical and  empirical  treatments  of  air  pollution,  air  concentration  is  found 
to  be  inversely  proportional  to  wind  speed.   Consequently,  the  air  quality 
downwind  of  a  freeway  must  be  closely  related  to  the  effective  height  of 
emission  and  elevated  highway  configuration  would  be  expected  to  result   in 
significantly  improved  downwind  air  quality  compared  to  at-grade  highways. 
Again,  we  must  qualify  this  statement  with  the  fact  that  highway  measurements 
adequate  to  verify  this  prediction  are  only  now  being  made. 

I  have  made  an  order  of  magnitude  calculation  of  the  wind  dilution  effect 
on  local  air  quality  to  be  gained  from  elevating  freeways.   Assuming  a 


logarithmic  wind  profile  and  all  other  factors  constant,  I  calculated  the 
effect  of  raising  the  highway  to  height  z  on  the  air  concentration  of  pollution, 
C,  at  the  same  height  at  some  arbitrary  distance  downwind.   I  obtained  the  re- 
lationship 

100  3C   -100 


C   3z   zlnz/z 
o 


.(3) 


If  we  take  z  =  0.1  m,  a  value  typical  of  suburban-land  use,  this  formula  gives 
the  following  results.   (See  Table  1-1) 

As  can  be  seen,  large  improvements  in  air  quality  would  be  expected  down- 
wind for  the  first  one  to  two  meters  of  elevation  but  with  decreasing  effect 
after  that.   It  should  be  emphasized  that  this  is  only  a  rough  order-of -mag- 
nitude estimate  and  neglects  other  effects  which  may  be  more  important  in  some 
cases,  but  would  be  expected  to  work  in  the  same  sense. 

A  second  qualifying  remark  should  be  made  at  this  point.   The  logarithmic 
wind  law  is  strictly  valid  only  for  conditions  in  which  the  atmosphere  is  well 
mixed.   Under  other  conditions,  corrections  must  be  made  to  this  relationship 
(Morgan,  Pruitt  and  Lourence,  1971)  which  are,  however,  not  large  in  magnitude 
and  are  well  known  providing  a  parameter  called  the  Richardson  Number,  which 
we  will  discuss  in  the  following  section,  is  measured  or  estimated. 

At  the  same  time  that  highway  pollution  is  being  transported  by  the  wind, 
it  is  being  diffused  by  atmospheric  turbulence.   It  is  useful  to  consider  the 
energetics  of  turbulence  at  this  point.   We  may  define  the  kinetic  energy  per 
unit  mass,  e,  of  the  field  of  turbulence  at  a  point  as 


e  =  1/2  (u)2  ,(4) 

where  u  is  the  turbulent  component  of  the  wind  speed  near  the  ground  such  that 

U  =  U  +  u  ,(5) 

where  U  is  the  average  wind  speed.   In  both  equation  (4)  and  (5) ,  the  overbar 

signifies  a  spatial  or  a  sufficiently  long  time  average.   Since  we  have  chosen 

to  consider  turbulence  as  energy,  it  is  instructive  to  look  at  the  sources  and 
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sinks  of  turbulence  energy  neat  the  ground.   This  problem  was  first  treated  by 
Richardson  (1920)  who  was  particularly  interested  in  the  special  case  when 
turbulence  vanishes  in  the  atmosphere,  a  situation  of  critical  significance  to 
the  air  pollution  problem. 

The  rate  at  which  turbulence  energy  is  generated  mechanically  by  wind 
shear  associated  with  surface  roughness  is  given  by  (Lumley  and  Panofsky,  1964) 

Ojf-I   (f)2  ,(6) 

*  3z    m   3z 

where  K  is  the  turbulent  diffusivity  for  momentum  in  the  atmosphere  and  arises 
m 

from  the  definition  of  U^.   The  rate  at  which  turbulence  is  suppressed  by 
stable  temperature  gradients  near  the  ground  is  given  by 

T  pc    T  \   3z   T  \   3z  'K,J 

P 

where  g  is  the  acceleration  of  gravity,  T  is  air  temperature  in  degrees  Kelvin, 
c   is  the  specific  heat  capacity  per  unit  mass,  H  is  the  flux  of  sensible  heat 
carried  upward  by  turbulence  (negative  for  stable  conditions  in  which  heat  is 
diffused  downwards) ,  K,  is  the  turbulent  diffusivity  for  heat  and  0  is  the  pt- 
tential  temperature,  which  in  most  practical  cases  may  be  replaced  by  temper- 
ature.  It  is  apparent  that  the  sign  of  this  term  changes  if  the  sensible  heat 
flux  is  upward,  which  corresponds  to  the  normal  condition  during  the  day  when 
heat  is  being  transferred  upwards  by  convection  from  the  warm  ground.   In  this 
case,  the  temperature  stratification  is  unstable  and  turbulence  is  being  gen- 
erated by  buoyant  forces  rather  than  being  suppressed.   For  stable  stratifi- 
cation, Richardson  reasoned  that  turbulence  would  disappear  when 

3T   „   ,3u\2 


T  %   3z    m  (3z; 


or 


T   n  3z    .  /oN 

m  3z 


that  is,  when  the  turbulence  sink  exceeds  the  source.   The  Richardson  number, 
Ri,  has  come  to  be  defined  by 

Ri  =  Tlfr-  '(9) 

where  potential  temperature  must  be  used  in  the  numerator  if  large  height  incre- 
ments are  used  in  evaluating  the  temperature  gradient.   Experimentally,  it  has 
been  found  that  turbulence  vanishes  at  Richardson  numbers  around  0.3  (Townsend, 
1957) .   That  fact  that  the  critical  Richardson  number  is  not  unity  reflects 
experimental  (Morgan,  Pruitt  and  Lourence,  1971)  observations  that  the  ratio 
K,  /K  is  considerably  less  than  one  under  stable  conditions. 

In  the  years  since  Richardson's  pioneering  work,  the  parameter  which  bears 
his  name  has  assumed  an  overwhelming  significance  in  micrometeorology .   It  has 
been  shown  experimentally  (Deacon,  1949)  and  theoretically  (Batchelor,  1953) 
that  the  Richardson  number  is  the  most  important  single  parameter  in  all  micro- 
meteorological  processes.   That  is  to  say,  Ri  is  more  important  than,  say, 
wind  speed  or  the  temperature  stratification,  considered  alone.   In  particular, 
the  Richardson  number  has  been  found  to  be  enormously  useful  in  organizing 
data  into  simple  and  understandable  patterns.   Figure  1-2  shows  some  typical 
micrometeorological  data  organized  into  useful  functional  relationships  by 
the  use  of  the  Richardson  number. 

The  Richardson  number  has  the  disadvantage  that  its  numerical  value  varies 
with  height  so  that  measurements  of  this  parameter  must  be  made  at  a  standard 
height  to  be  comparable.   Another  equivalent  parameter,  more  useful  in  some 
ways,  may  be  obtained  by  forming  the  ratio  of  the  two  energy  source/sink  terms 
as  before  but  now  using  the  most  fundamental  definition  on  the  left  of 
equations  (6)  and  (7) , 
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where  we  have  used  the  vertical  derivative  of  the  logarithmic  wind  law  in  this 

transformation.   The  quantity  in  brackets,  which  has  the  dimension  of  length 

and  is  approximately  constant  with  height,  is  called  the  Monin-Obukhov  length, 

L,  where  » 

"  P^  T  U* 
L=— k^ '(11) 

Under  stable  conditions  (downward,  negative  heat  flux,  H) ,  the  Monin-Obukhov 
length  has  the  interpretation  of  being  roughly  the  height  at  which  turbulence 
is  suppressed.   Under  stable  conditions,  then,  pollution  emitted  near  the 
ground  would  be  expected  to  diffuse  to  height  L.   Under  very  stable  conditions, 
L  becomes  small  and  air  pollution  concentrations  high.   Under  unstable  con- 
ditions, L  is  negative  and  has  the  physical  interpretation  as  being  the  height 
at  which  convectively-produced  turbulence  energy  compares  with  mechanically 
produced  energy.   Above  L,  convection  predominates  and,  under  these  conditions, 
turbulence  levels  are  higher,  atmospheric  diffusion  more  efficient,  and  air  con- 
centrations of  pollution  are  less. 

We  shall  discuss  the  mathematics  of  diffusion  in  the  next  section.   In 
qualitative  terms,  it  is  permissible  to  visualize  atmospheric  dispersion  as  a 
process  similar  to  molecular  diffusion  in  a  solid  or  liquid  body.   The  random 
motion  of  the  molecules  acts  on  superimposed  gradients  to  transport  properties, 
such  as  heat,  "down  the  gradient"  from  regions  of  high  concentration  to  low. 
In  the  analogy  with  turbulent  motion  in  fluids,  turbulence  is  thought  of  as 
essentially  random  motion  which  acts  to  mix  or  equalize  distribution  of  fluid 
properties  and,  thus,  the  net  result  is  a  transport  process  in  which  one  re- 
gion of  the  fluid  gains  at  the  expense  of  another.   Formally  we  can  define 
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dif fusivity ,  K  ,  of  a  fluid  property,  s,  expressed  in  units  per  unit  mass  of 


fluid  by 


F   =-pK  f  ,(12) 

sx       s  8x 


where  F   is  the  flux   (transport  per  unit  area  and  time)  of  s  in  the  x  dir- 
ection resulting  from  the  diffusion  process,  and  c  is  the  concentration  per 
unit  mass  of  air.   The  dif fusivity  so  defined  is  formally  equivalent  to  mole- 
cular dif fusivity  and,  indeed,  may  be  regarded  as  the  sum  of  the  molecular  and 
turbulent  processes.   The  size  of  the  turbulent  component  of  the  dif fusivity 
is,  however,  many  orders  of  magnitude  larger  than  the  molecular  term  in 
natural  fluid  systems.   The  conclusion  to  be  drawn  is  then  that  turbulence  is 
vastly  more  efficient  at  transporting  fluid  mass  and  whatever  properties  as- 
sociated, including  pollution,  than  is  molecular  diffusion.   That  is  the 
central  significance  of  turbulence  and  also,  as  far  as  we  know,  the  fundamental 
reason  for  its  existence. 

The  molecular  analogy  is  useful  only  in  a  very  general  sence.   Unlike 
molecular  properties,  such  as  conductivity  or  viscosity  which  can  be  tabulated 
as  physical  constants  for  many  systems,  turbulent  dif fusivity  is  a  complex 
function  of  the  state  of  the  flow.   It  is,  for  instance,  a  function  of  the 
Richardson  number.   The  property  of  turbulence  which  makes  it  particularly 
difficult  to  deal  with  both  practically  and  theoretically  is  that  the  dif- 
fusivity  cannot,  in  principle,  be  treated  as  constant  even  in  situations  where 
the  measured  Richardson  number  is  constant.   The  reason  is  that  dif fusivity  is 
observed  to  be  a  function  of  the  scale  of  fluid  motion.   If  one  follows  a  puff 
of  smoke  emitted  from  a  source  near  the  ground,  it  is  easy  to  observe  a  rapid 
increase  of  size  of  the  puff  as  it  entrains  (mixes  with)  fluid  from  the  en- 
vironment.  This  process  is  a  function  of  the  energy  which  is  contained  in 
scales  of  motion  (or  "eddies")  which  are  smaller  than  the  puff.   It  is  a 

fundamental  property  of  the  atmosphere,  and  all  turbulent  flow  systems,  that 
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if  measurements  are  made  of  the  amount  of  kinetic  energy  available  in  the 
various  scales  of  motion  present  in  the  turbulence,  it  is  invariably  found 
that  more  energy  exists  at  large  scales  than  at  small  scales.   Consequently, 
as  our  smoke  puff  grows,  it  is  subjected  to  more  and  more  energetic  diffusion 
by  the  scales  of  turbulent  motion  smaller  than  its  current  dimensions.   The 
larger  it  gets,  the  faster  it  diffuses.   L.  F.  Richardson  (1926)  was  the  first 
to  recognize  this  remarkable  phenomenon,  and  he  proposed  that  this  variation 

could  be  well  represented  by 

4/3 
K  -  constant  x  1  ,(13) 

s 

where  1  is  a  representative  size  of  the  "eddy"  or  identifiable  diffusing  sub- 
stance.  "Richardson's  law"  has  stood  the  test  of  time  and  now  constitutes  a 
primary  objective  of  any  new  turbulence  theory. 

There  are  several  consequences  of  these  properties  of  turbulence  which 
are  of  significance  to  the  practical  objective  of  modeling  the  dispersion  of 
pollution  in  the  atmosphere.   If  one  clooses  to  model  the  dispersion  of  ef- 
fluent by  following  individual  "puffs"  or  identifiable  portions  of  the  polluted 
air,  it  is  not  possible  to  consider  the  diffusivity  as  constant.   Richardson's 
law  must  be  taken  into  account.   Second,  if  one  considers  the  distribution  of 
turbulence  energy  with  height  above  the  surface  of  the  earth,  it  is  reasonable 
to  expect  that  as  you  gain  altitude,  larger  and  larger  scales  of  motion  will 
exist  because  there  is  more  "room".   This  has  been  found  to  be  true  observa- 
tionally  (Myrup ,  1967)  along  with  the  logical  corollary  that  diffusivity  in- 
creases strongly  with  height.   Hence,  if  the  Richardson  number  at  a  given 
height  is  constant  in  time,  it  is  permissible  to  consider  the  diffusivity  as 
constant  at  that  height.   However,  if  the  process  extends  over  a  considerable 
range  in  height,  as  diffusion  does,  the  diffusivity  is  not  a  constant  but  a 
function  of  height.   Thus,  the  common  assumption  in  air  pollution  meteorology 

that  the  diffusivity  is  a  constant  quantity  is  an  approximation  whose  limitations 
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should  be  explicitly  recognized. 

At  this  point,  I  wish  to  point  out  that  the  dispersion  of  particulates, 
large  enough  to  fall  out  near  their  source,  has  certain  unique  properties. 
Since  particulates  emitted  from  any  source  have  a  spectrum  of  sizes,  the  fall- 
out is  differential;  big  particles  reach  the  ground  first.   This  effect  alone, 
in  the  absence  of  any  turbulence  or  wind  shear,  is  enough  to  produce  a  wide 
dispersion  of  the  material  which  falls  out  and  is  deposited  at  the  surface  of 
the  earth.   This  effect  is  illustrated  in  Table  1-2  in  which  we  list  the  ter- 
minal velocity  for  various  sized  particles  and  a  single  order-of-magnitude 
calculation  of  the  distance  downwind  at  which  a  particle  of  the  indicated  size 
would  be  deposited  if  the  average  wind  speed  was  2  ms  "  and  the  release  height 

was  2  meters.   As  can  be  seen,  the  dispersion  due  to  differential  fallout 
is  considerable.   This  is  also  a  stability-dependent  process,  since  the  shape 
of  the  wind  profile  (and  hence  the  average  value  of  the  wind)  and  the  turbulent 
diffusivity  between  the  release  height  and  the  surface  are  a  function  of  the 
Richardson  number. 

The  magnitude  of  the  differential  fallout  effect  is  heavily  dependent  on 
the  shape  of  the  particulate  size  spectrum.   To  my  knowledge,  insufficient 
measurements  have  been  made  of  the  size  distribution  of  particulates  emitted 
by  motor  vehicles  to  allow  modeling  of  this  process  with  any  degree  of  con- 
fidence. 

The  horizontal  dimensions  of  cities  are  such  that  many  properties  of  the 
atmosphere  above  them  belong  to  those  of  the  meteorological  meso-scale.   The 
term  "air  pollution  meteorology"  usually  refers  to  the  phenomena  of  the  meso- 
scale  and,  indeed,  most  air  pollution  studies  have  been  concerned  with  this 
scale.   The  reason  for  making  this  distinction  is  that  many  of  the  fundamental 
properties  of  this  scale  of  motion  are  different  from  those  of  the  micro-scale. 
Table  1-3  is  a  schematic  summary  of  the  contrasting  properties  of  the  two  scales 
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TABLE  1-2. 

Particle  radius (y)    Terminal  velocity(cms  )     Deposition  point 
downwind (m) 

5                       .3  133.3 

10                      1.3  30.8 

50                     32.0  12.5 

100                    136.0  2.94 

400                   340.0  1.18 
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TABLE   1-3 


Microscale 

Mesoscale 

Phenomena 

Pollutant 

Primary  gaseous  and  particulate  pollutant 

i 

Products  of  photochemical 
reactions 

Space  scale 

1  -  100  m 
(highway  corridor) 

0.1  to  100  Km 
(city,  air  basin) 

Time  scale 

1-60  minutes 

1  hr.  to  1  day 

Primary  Transport 

Prevailing  (geostrophic) 

wind 

processes 

Turbulence 

Diurnal  wind  systems  (sea 
breezes,  etc.) 

Source  of  transport 
energy 

Surface  roughness, 
wind  shear, 
convection 

Horizontal  temperature 
contrasts,  topography 

Primary  parameters 

Wind  speed,  Surface  roue 

hness 

Richardson  Number 

Latitude,  stability 
(inversion  height) 

Typical  model 

Line-source 
Goussian 

Box 
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from  which  their  differing  significances  to  pollution  derive.   Meso-scale  wind 
systems  arise  because  of  horizontal  gradients  in  the  temperature  of  the  earth's 
surface.   Sea  and  lake  breezes  are  familiar  examples  of  this  phenomenon.   These 
wind  systems  have  the  properties  that  the  strongest  wind  speeds  are  usually 
near  the  surface  of  the  earth,  and  wind  direction  often  exactly  reverses  with 
height.   Figures  1-3  and  1-4  show  examples  of  pollution  transport  in  the  Los 
Angeles  basin  by  the  sea  breeze,  and  Figure  1-5  illustrates  the  typical  re- 
versal of  wind  direction  with  height.   Since  meso-scale  wind  systems  are 
strongly  associated  with  the  24-hour  period  of  the  sun,  the  wind  flow  develops 
slowly  enough  that  the  component  of  the  earth's  rotation  about  the  local  ver- 
tical (Coriolis  parameter),  which  is  a  function  of  latitude,  becomes  an  in- 
fluence.  Measurements  of  the  Richardson  number  are  difficult  to  make  in  the 
free  atmosphere.   Consequently,  atmospheric  stability  is  usually  specified  by 
means  of  the  lapse  rate  alone.   Elevated  inversion  (temperature  increase  with 
height)  layers  are  usually  present  in  diurnal  wind  systems  and  the  height  of 
the  base  of  the  inversion  is  an  important  parameter  for  meso-scale  air  pol- 
lution modeling.   Within  inversion  layers,  the  Richardson  number  often  be- 
comes large  enough  that  all  turbulence  ceases,  and  hence  inversions  are  ef- 
fective barriers  to  atmospheric  dispersion.   Figure  1-6  is  a  typical  temper- 
ature sounding  in  the  Los  Angeles  basin  and  shows  the  famous  inversion  layer 
which  is  usually  present  over  that  area  during  the  summer  and  fall  months. 

In  qualitative  summary,  we  may  say  that  the  most  important  physical  pro- 
cesses and  parameters  which  determine  the  concentration  of  pollution  in  the 
air  which  we  breathe,  are  the  following  (not  necessarily  in  order  of  im- 
portance) :   rate  at  which  pollution  is  emitted  into  the  atmosphere("source 
strength"),  distance  downwind  from  the  source  of  pollution,  wind  speed,  at- 
mospheric stability  near  the  ground  (Richardson  number) ,  height  of  elevated 
inversion  layers  and  topography. 
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M.  NEIBURGER 


Figures  1-3  and  1-4.  The  Los  Angeles  sea  breeze. 
(Neiburger  et  al,  1956) . 
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NORTH-SOUTH  WIND  SPEED 
SACRAMENTO,  CALIFORNIA 
AVERAGE  OF   3  OBSERVATION 
PERIODS       JULY-AUGUST  1971 


NORTH -SOUTH  WIND    COMPONENT    (ids'1) 
Figure  1-5.     Typical  meso-scale  wind-speed  profile. 
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Figure  1-6.  Typical  sounding  of  tenperature  and 
humidity  over  Los  Angeles. 
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Mathematical  Modeling  of  Air  Quality 

In  this  section  we  shall  discuss  the  application  of  mathematical  and 

numerical  techniques  to  the  problem  of  estimating  air  quality  given  information 

on  pollution  emission.   With  reasonable  generality,  we  can  state  the  governing 

equation  for  the  concentration,  C,  of  some  pollutant  as 

9C  9C  9C  9C  9  ,_  3Cv  9  ,v  3C.  9  ,_  9C. 
—  =  -u  —  -v  —  -w  —  +  —  (K  — )  +  —  (K  — )  +  — -  (K  — ) 
9t      9x    9y    9z    9x   x  9x    9y   y  9y    9z    z  9z 

+R(t)C  +S(x,y,z,t)  ,(14) 

where  u  and  v   are  the  horizontal  and  w  the  vertical  components  of  the  velocity 

of  the  air  in  the  x,y,  and  z  directions,  K  ,  K  ,  and  K  are  turbulent  diffusiv- 

J '  '  x   y       z 

ities  associated  with  turbulent  transport  in  the  three  coordinate  directions, 
R(t)  is  a  chemical  reaction  function  for  the  case  in  which  the  pollutant  is 
chemically  active,  and  S(x,y,z,t)  is  a  completely  artibrary  source  function 
which  expresses  the  rate  at  which  pollutant  is  emitted  into  a  unit  volume. 
Many  different  approaches  can  be  adopted  in  obtaining  analytic  or  numerical 
solutions  to  this  equation.   We  shall  concentrate  on  the  most  common  solutions 
and  methods  which  are  currently  in  operational  use  in  estimating  air  quality. 

As  it  stands,  equation  (14)  is  formidable  and  solutions  are  only  possible 
if  it  is  greatly  simplified.   For  instance,  if  the  mean  wind  is  in  the  x  dir- 
ection, the  vertical  velocity  is  neglected,  diffusion  in  the  x  direction  is 
neglected  relative  to  wind  transport,  the  pollutant  is  non-reacting,  and  the 
dif fusivities  are  taken  as  constants,  the  equation  (10)  reduces  to 

f  =-Uf  +Ky^f+Kr^|+S(x.,...t)  .(15) 

9y       9z 

For  this  case,  application  of  boundary  conditions  appropriate  for  certain 

source  configuration  allows  explicit  analytic  solutions  to  be  found.   Equation 

(15)  is,  in  fact,  essentially  the  molecular  or  Fickian  diffusion  equation  and 

numerous  solutions  are  available  in  any  text  on  heat  conduction  or  molecular 

diffusion  (Crank,  1956).   For  instance,  one  of  the  most  commonly  used  solutions, 
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which  is  valid  for  a  point  source  of  pollution  emitted  at  height  h  is 


C  = 


2tt  a  a  U 
y  z 


Exp 


2  V 


.Exp 


1  ,z-h.2 

2  Caz  ; 


Exp 


i  /-z±h\2 
2  laz  ; 


(16) 


where  Q  is  the  rate  of  release  of  the  pollutant,  U  is  the  average  wind  speed, 

and  a  and  a     are  the  standard  deviations  of  the  concentration  from  the  plume 
y      z  r 

axis  as  functions  of  distance  downwind.   This  formula  is  formally  analogous  to 
normal  or  Gaussian  bivariate  probability  distribution,  and  the  term  "Gaussian" 
is  usually  applied  to  models  which  use  equations  of  this  sort.   Turner  (1963) 
used  this  solution  in  an  air  pollution  model  for  computing  the  sulfur  dioxide 
concentration  in  Nashville,  Tennessee. 

The  key  to  success  for  this  or  any  other  Gaussian  formula  is  selection  of 
the  appropriate  a's  as  functions  of  distance  downwind.   These  terms  are,  in 
fact,  functions  of  the  diffusivity  and  hence  functions  of  surface  roughness 
and  atmospheric  stability,  as  measured  by  the  Richardson  number.   Turner 
(1969)  has  published  a  practical  system  for  making  air  quality  estimates  based 
on  procedures  developed  by  Pasquill  (1961)  and  Gifford  (1961) .   In  this  scheme 
stability  is  estimated  from  meteorological  conditions,  such  as  wind  speed, 
cloud  cover,  and  intensity  of  solar  radiation,  in  terms  of  six  stability 
classes,  A(strongly  unstable)  to  D  (neutral)  to  F(strongly  stable).   The  dis- 
persion parameters,  a  and  a  are  then  taken  from  graphs  which  specify  down- 

y     z 

wind  variation  as  a  function  of  stability.   Figure  1-7  shows  the  set  of  curves 

for  K  . 
z 

Equation  (16)  includes  the  assumption  that  the  surface  acts  like  a  per- 
fect reflector,  i.e.,  no  deposition  or  absorption  takes  place.   Clearly  this 
assumption  is  appropriate  only  for  chemically  inactive  gases.   The  presence 
of  elevated  inversion  layers  can  be  handled  by  assuming  that  when  the  plume 
dimensions  have  become  comparable  to  the  height  of  the  base  of  the  inversion, 
H,  this  becomes  a  lid  preventing  further  vertical  growth  and  thereafter  mixing 
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Figure  1-7.  Vertical  dispersion  coefficient  as  a  function  of 
downwind  distance  from  the  source. 
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in  the  vertical  is  complete.   If  this  distance  downwind  is  assumed  to  be  twice 
that  at  which  one-tenth  of  the  plume  has  penetrated   the  inversion,  then  the 
formula  for  air  concentration  becomes 

,(17) 


C  =  2 _  Exp 

/2tT  a  H  U 

y 


2  Ko  ' 

y 


where  the  concentration  is  now  independent  of  the  vertical  direction  and  this 
relation  is  supposed  to  be  valid  when  a  >  0.94H. 

The  Gaussian  approach  is  easily  extended  to  other  source  configurations. 
For  a  continuously  emitting  infinite  line  source,  the  solution  at  ground  level 


is 


C £2 Exp 

v2tT  a  U  sin* 
z 


2  Ko  ' 
z 


,(18) 


where  <f>  is  the  angle  between  the  wind  direction  at  the  line  (zero  at  right 
angles)  . 

Equations  (16)  and  (18)  are  the  working  formulae  for  the  great  majority 
of  practical  schemes  for  estimating  air  concentration  from  air  pollution 
source  inventory  information.   In  order  to  calculate  the  air  concentration  at 
a  specific  point  in  an  urban  area,  the  city  is  divided  up  into  a  set  of  equal 
areas  defined  by  a  grid  and  an  appropriate  source  strength  assigned  to  each. 
For  instance,  in  his  Nashville  study,  Turner  obtained  S0„  source  strengths 
for  each  square  mile  in  a  17  x  16  mile  rectangle.   Downwind  concentrations  re- 
sulting from  S0„  emitted  within  a  given  square  mile  were  calculated  from  the 
point  source  formula,  assuming  that  the  area  emission  was  concentrated  at  a 
point.   This  calculation  was  repeated  for  each  of  272  square  mile  sources  and 
the  results  summed  to  give  the  net  result.   Observed  wind  speeds  and  a's  es- 
timated from  observed  meteorological  conditions  were  used.   Comparison  between 
observed  and  computed  24-hour  averages  showed  that  58%  of  all  computed  values 
were  within  -1  pphm,  a  result  which  was  taken  as  demonstrating  the  feasibility 

of  the  approach. 
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Several  other  investigators  have  made  calculations  similar  to  those  of 
Turner.   Koogler,  et  al.  (1967)  used  Turner's  equation  with  the  addition  of  an 
exponential  decay  to  account  for  S0„  removal  processes  for  Jacksonville,  Fla. 
Koogler  obtained  95%  of  his  calculated  SO-  concentrations  within  -1  pphm  of 
observed  values.   Hilst  (1967)  used  essentially  the  same  approach  for  the  state 
of  Connecticut  using  an  improved  method  of  incorporating  observed  winds.   5,600 
square  area  sources  were  used,  each  5,000  feet  on  a  side,  and  important  in- 
dividual sources,  such  as  large  power  plants,  were  treated  individually. 

The  most  sophisticated  air  pollution  model  now  in  existence  is,  to  our 
knowledge,  that  which  is  being  developed  by  Neiburger  and  Lamb  (1968).   This 
model  includes  improved  solutions  to  the  diffusion  equation  for  point,  finite 
line,  and  area  sources,  absorption  of  pollutants  at  the  ground,  simple  chem- 
ical reactions  whose  rates  may  be  given  as  arbitrary  functions  of  time,  and 
horizontally-  and  time-varying  winds.   In  this  model,  the  diffusion  coef- 
ficients are  constants  and  wind  does  not  vary  with  height,  and  vertical  motion 
is  neglected.   Preliminary  calculations  with  this  model,  using  observed  winds 
and  traffic  data,  were  made  for  carbon  monoxide.   The  results,  while  considered 
encouraging,  revealed  a  major  source  of  error  associated  with  the  neglect  of 
vertical  motion.   Unrealistically  large  values  of  carbon  monoxide  concentration 
were  calculated  to  occur  in  regions  of  horizontal  convergence  of  the  wind  in- 
dicating that  removal  of  pollutant  by  vertical  motion  is  an  important  process 
which  cannot  be  neglected.   This  model  is  in  an  active  state  of  development 
and  important  and  high  quality  results  are  anticipated  in  the  near  future. 

Since  we  can  illustrate  most  of  the  points  we  wish  to  make  by  reference 
to  the  models  we  have  already  discussed,  we  shall  limit  our  presentation  of 
the  remaining  major  approaches  to  brief  outlines.   In  addition  to  the  Gaussian 
models  and  Neiburger1 s  and  Lamb's  improved  version  of  this  technique,  there 
are  four  other  major  approaches  to  the  problem  of  calculating  air  quality 
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from  measurable  parameters.   One  is  based  on  the  statistical  theory  of  tur- 
bulence, originally  introduced  by  G.  I.  Taylor  (1921).   In  this  approach, 
emphasis  is  placed  on  the  statistical  properties  of  velocity  fluctuations  of 
ensembles  of  passively  floating  particles  or  marked  fluid.   The  statistical 
theory  deals  directly  with  the  spread  of  a  group  of  particles  with  respect  to 
a  frame  of  reference  moving  with  the  fluid  (Lagrangian  frame)  which  is  the  most 
fundamental  manner  in  which  to   treat   turbulent  diffusion.   Since  measurements 
are  almost  always  made  with  respect  to  a  fixed  (Eulerian)  frame  of  reference, 
a  central  problem  in  applying  the  statistical  theory  to  practical  dispersion 
problems  is  to  establish  the  correspondence  between  Lagrangian  and  Eulerian 
turbulence  parameters.   The  statistical  theory  may  be  the  most  fundamental  but 
is  probably  the  most  remote  from  immediate  practical  application  at  this  time. 

Another  major  approach  which  is  only  now  being  developed  and  hopefully 
represents  a  useful  compromise  between  the  unashamed  empiricism  of  the  Gaussian 
models  and  more  satisfying  basic  theory  is  the  hydrodynamics  approach  (Knox, 
et  al. ,  1971).   In  this  technique,  the  fundamental  equations  of  meteorological 
dynamics  and  thermodynamics  are  numerically  solved  on  a  three-dimensional  grid 
using  approximate  topography  as  a  boundary  condition  to  provide  the  mesoscale 
flow  and  stability.   Smaller  scale  flow  features,  or  "turbulence  ",  are  handled 
by  means  of  diffusivity  parameters,  such  as  we  have  been  discussing.   A  less 
satisfying  variant  of  this  approach  would  be  to  use  observed  three-dimensional 
winds  to  specify  the  mesoscale  flow.   This  alternative  appears  less  satisfying 
to  us  because  of  the  extraordinary  measurement  problem. 

Three-dimensional  numerical  integrations  are  so  demanding  of  computer 
storage  and  computational  speed  that  only  relatively  crude  models  can  be  solved 
at  this  time.   It  remains  to  be  seen  whether  computer  capabilities  develop 
fast  enough  for  this  approach  to  be  of  practical  utility  in  the  near  future. 

The  next  of  the  major  approaches  to  our  problem  is  the  similarity  theory, 
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introduced  by  Monin  (1959)  and  developed  further  by  Batchelor  (1959,  1964), 
Gifford  (1962)  and  Pasquill  (1966).   In  this  technique,  a  semi-empirical  frame- 
work based  on  dimensional  analysis  is  adopted  to  organize  pollutant  data  in 
terms  compatible  to  micrometeorological  practice.   We  present  some  results  in 
the  next  section  obtained  with  what  is  essentially  a  similarity  approach,  so 
we  will  not  develop  this  approach  further  at  this  time. 

The  last  approach  to  be  discussed  here  is  the  so-called  box  model,  which 
is  almost  trivially  simple  but  nevertheless  useful  for  mesoscale  applications. 
In  situations  where  a  well-defined  inversion  acts  as  an  effective  lid  on  ver- 
tical dispersion  and  when  pollution  is  emitted  more  or  less  uniformally  over  a 
large  urban  area,  it  may  often  be  acceptable  to  calculate  concentration  from 

C=^  ,(19) 

UH 

where  U  is  the  average  wind  speed  between  the  surface  and  the  height  of  the 
base  of  the  inversion,  k  is  a  constant,  and  Q  is  an  emission  rate  per  unit  area. 
Hanna  (1971)  has  shown  that  simple  formulae  like  equation  (19)  are  often  as 
accurate  as  much  more  sophisticated  approaches,  as  near  as  can  be  told  from 
the  available  data. 

The  modeling  approaches  we  have  been  discussing  are  meant  to  be  used  in  a 
practical,  engineering  sence.   They  must  ultimately  be  judged  in  this  spirit. 
Before  attempting  an  overall  evaluation  of  the  state  of  the  art  of  highway  air 
quality  modeling,  it  is  instructive  to  look  at  an  operational  attempt  to  use 
some  of  the  techniques  we  have  presented. 

In  the  last  two  years  the  California  Division  of  Highways  has  been  con- 
ducting an  extensive  environmental  impact  program  to  evaluate  the  effect  on  air 
quality  of  current  highways  in  the  state  and  to  make  projections  of  the  effect 
of  future  planning  decisions  regarding  highway  construction.   At  present,  this 
program  is  limited  to  an  evaluation  of  the  probable  air  quality  impact  of  the 
decision  to  build  or  not  to  build  highways  whose  routes  and  configurations  were 

planned  well  before  the  overriding  significance  of  environmental  impact  became 
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clear.   There  were  four  major  objectives  for  the  program.   One  was  to  generate 
reliable  data  on  the  current  state  of  air  quality  near  highways.   A  second  ob- 
jective was  to  develop  and  verify  practical  highway  air  quality  models.   The 
third  was  to  develop  models  for  assessing  highway  impact  on  the  air  quality  of 
entire  mesoscale  air  basins.   A  fourth  was  to  use  verified  models  to  estimate 
air  quality  in  twenty  years  time  for  the  built  and  not-built  cases  for  par- 
ticular highway  plans.   At  present,  several  contracts  have  been  awarded  to 
research  and  development  groups  in  the  private  sector  to  perform  these  impact 
studies.   I  shall  briefly  outline  the  results  of  the  first  of  the  studies  to 
be  completed,  which  was  an  air-quality  analysis  and  impact  study  of  proposed 
routes  92/238  near  Hayward ,  California.   The  Hayward  area  lies  on  the  coastal 
plain  to  the  east  of  the  southern  portion  of  the  San  Francisco  Bay.   The  im- 
mediate topography  is  uncomplicated;   on  a  larger  scale,  Hayward  lies  in  the 
basin  which  encloses  the  Bay.   The  firm  which  was  awarded  the  contract  col- 
lected all  available  air  quality  and  meteorological  data  for  the  area,  es- 
tablished additional  meteorological  and  air  pollution  monitoring  stations, 
developed  a  Gaussian  transport  and  diffusion  model,  verified  the  model  with 
data  taken  during  a  special  concentrated  observation  period,  and  performed 
calculations  with  the  model  to  estimate  air  quality  in  the  region  for  the  year 
1990  with  and  without  the  proposed  freeways.   On  the  whole,  there  is  no  doubt 
that  this  was  a  highly  professional  piece  of  applied  science  and  is  well  re- 
presentative of  the  current  state  of  our  measurement  techniques  and  modeling 
capability.   For  these  reasons,  whatever  criticism  we  make  of  the  techniques 
employed  or  results  obtained  in  this  study  are  to  be  taken  as  applying  to  the 
field  as  a  whole. 

Figure  1-8  shows  an  example  of  the  basic  data  which  was  collected  in  the 
Hayward  area  by  the  contractor.   The  pollutant  is  carbon  monoxide  averaged  over 
24  hours.   The  numbers  attached  to  the  data  points  refer  to  the  various  sites 
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Figure  1-8.  Time  variation  of  carbon  monoxide  concentration 
at  various  sites  in  the  vicinity  of  Hayward,  California  in 
1971. 
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used  in  the  study.   The  extreme  day-to-day  variability  is  typical  of  pollution 
data.   The  horizontal  dashed  line  is  the  ambient  air  quality  standard  for 
carbon  monoxide  (  12  hours  at  10  ppm) .   Figure  1-9  shows  a  cumulative  frequency 
diagram  derived  from  the  1063  hourly  measurements  of  CO  made  during  the  two 
and  a  half  months  of  this  study.   The  diagram  shows,  for  instance,  that  hourly 
average  CO  concentrations  of  10  ppm  were  equalled  or  exceeded  18%  of  the  time 
during  this  study.   Figure  1-10  shows  the  non-cumulative  frequency  distribution 
for  the  same  data  but  now  broken  down  by  site.   The  variability  from  site  to 
side  and  the  apparent  bimodal  character  of  these  curves  is  striking.   In  ad- 
dition, it  could  be  said  that  these  curves  do  not  have  the  usual  smooth  ap- 
pearance of  well-defined  and  stable  statistics.   The  question  must  be  asked  as 
to  whether  the  conditions  of  this  study  allowed  time  for  sufficient  data  to  be 
gathered  to  allow  reliable  statistical  generalizations  to  be  made.   Figure  1-11 
shows  a  similar  analysis  of  weekday  versus  weekend  data  for  all  sites.   The 
same  comments  apply. 

The  investigators  state,  although  the  analysis  is  not  presented  in  suf- 
ficient detail  to  allow  evaluation,  that  a  significant  correlation  was  found 
between  carbon  monoxide  and  the  product  UH,  where  these  parameters  have  the 
same  meaning  as  in  the  box- model  approach  discussed  above,  so  that  the  result 
amounts  to  be  partial  validation  of  that  approach. 

Figure  1-12  shows  a  summary  of  a  concentrated  microscale  study  in  the  im- 
mediate vicinity  of  an  existing  freeway  section  which  ran  approximately  North- 
South.   All  data  are  shown  from  a  variety  of  meteorological  and  traffic  situa- 
tions and  it  is  difficult  to  see  a  clear  pattern  such  as  simple  models  such  as 
the  line-source  Gaussian  would  lead  you  to  expect.   Figure  1-13  shows  a  typical 
individual  case  for  downwind  data  and  Figure  1-14  for  upwind  data.   While  these 
displays  appear  to  be  more  reasonable  than  the  previous  scatter  diagram,  serious 
methodological  questions  arise  from  consideration  of  these  data.   The  measurements 
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monoxide  concentration  at  Hayward,  California. 
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Figure  1-10.  Frequency  distribution  of  CO  concentration  for 
various  sites  at  Hayward,  California. 
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Figure  1-11.  Frequency  distribution  of  00  concentration  on 
weekdays  and  weekends  for  Hayward,  California. 
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Figure  1-12.     CO  concentrations  in  the  inmsdiate  vicinity 
of  the  Nimitz  Freeway  near  Hayward,  California.     All 
measurements  shown. 
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LINESITE     STUDY:     Km,  #s  uct.u  (Mo.,.)  0955-1412  (pst) 
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Figure  1-13.  Pollution  measurements  made  upwind  of  Nimitz 
Freeway  near  Hayward,  California,  Oct.  11,  1971. 
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LINESITE      STUDY:       Run  #7  Oct.    12   (Tues.)  1005-1635  (PST) 
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Figure  1-14.  Pollution  measurements  made  downwind  of 
Nimitz  Freeway  near  Hayward,  California.  Oct.  12,  1971. 


37 


were  made  by  moving  a  van  containing  air  pollution  sensing  apparatus  from  point 
to  point.   Thus,  the  data  are  not  simultaneous  and  contain  large  variations  of 
meteorological  parameters  from  point  to  point.   Such  a  procedure  would  not  be 
acceptable  in  any  professional  micrometeorological  study.   Measurements  must 
be  simultaneous  to  be  comparable.   I  have  no  doubt  the  contractor  did  the  best 
job  possible  with  the  time  and  resources  available.   Nevertheless,  it  must  be 
said  that  we  need  to  do  much  better  than  this. 

A  validation  analysis  was  made  of  the  highway  diffusion  model,  and  this  is 
shown  in  Figure  1-15.   The  objective  was  to  estimate  the  contribution  to  carbon 
monoxide  concentration  at  specific  locations  from  an  existing  freeway  alone. 
As  can  be  seen,  the  "background"  concentrations  are  larger  than  the  calculated 
freeway  contributions.   It  is  difficult  to  have  confidence  in  the  modeling  ap- 
proach on  the  basis  of  this  analysis  alone. 

An  additional  criticism  can  be  made  concerning  the  lack  of  analysis  of  the 
properties  of  the  model  itself.   At  a  minimum,  a  sensitivity  analysis  of  any 
operational  model  should  be  made  to  determine  the  relative  importance  of  the 
various  input  data  and  the  degree  of  precision  needed.   For  instance,  how  sen- 
sitive are  the  validation  calculations  discussed  above  to  uncertainties  in,  say, 
the  wind  speed  or  freeway  emission  rate?   Without  such  information,  it  is  dif- 
ficult to  interpret  a  validation  analysis. 

To  fulfill  the  final  portion  of  the  study,  the  contractor  used  the  de- 
veloped models  to  make  air-quality  projections  for  the  year  1990.   On  the  meso- 
scale  a  number  of  calculations  were  made  for  various  wind  directions  and  sta- 
bilities for  the  cases  in  which  the  proposed  routes  92/238  were  built  and  not- 
built.   In  general,  the  mesoscale  calculations  showed  little  difference  be- 
tween these  two  cases. 

On  the  microscale,  within  a  thousand  feet  of  the  highway  corridor,  re- 
latively large  differences  were  found  between  the  built  and  not-built  cases. 
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Figure  1-15.  Test  of  an  air  pollution  transport  and  diffusion 
model. 


Date NAWC  Model      Background      Computed Observed 

(A)  -  WSW  Flow 

9-30-71  0.2  1.4  1.6  1.5 

(1400-1700  PST) 

10-1-71  0.4  2.4  2.7  2.8 

(0800-1100  PST) 

10-21-71  0.3  2.8  3.1  3.8 

(1200-1600  PST) 

11-7-71  0.4  4.7  5.1  5.0 

(0900-1100  PST) 


(B)  -  WSW  Flow 

10-29-71  1.5  3.5  5.0  6.5 

(1400-1500  PST) 


(C)  -  Drainage  Flow 

10-3-71  1.9  1.2  3.1  4.6 

(0400-0700  PST) 

11-4-71 
(0400-0700  PST)  2.4  2.4  4.4  8.4 
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Figure  1-16  shows  typical  results  which  indicate  lower  upwind  concentrations 
of  carbon  monoxide  and  considerably  higher  downwind  concentrations  of  carbon 
monoxide  for  the  built  case  in  comparison  with  not-built.   In  view  of  the  un- 
certainties in  the  validation  program  discussed  above,  it  becomes  doubly  un- 
certain as  to  how  to  interpret  or  use  these  microscale  projections.   An  im- 
portant point  to  note  is  that  these  calculations  are  presented,  as  is  appropri- 
ate, in  the  form  of  air  concentrations  at  various  distances  from  the  freeway. 
The  validation  measurements,  however,  were  made  only  at  one  point,  and  we  have 
no  way  of  knowing  the  ability  of  the  model  to  predict  spatial  variation. 

It  is  difficult  to  assess  the  overall  state  of  the  field  of  air  pollution 
modeling.   The  field  is  moving  very  fast  and  many  complex  problems  are  involved 
and  must  be  considered.   In  particular,  the  measurement  phase  is  critically  im- 
portant and  must  be  discussed  in  connection  with  the  modeling  effort. 

First  of  all,  it  is  clear  the  physical  rigor  of  all  operational  and  re- 
search models  is  not  high.   The  treatment  of  diffusion  parameters,  for  in- 
stance, is  usually  highly  arbitrary  and  unconnected  with  what  we  know  of  the 
physics  and  energetics  of  turbulence.   In  particular,  diffusion  coefficients 
are  almost  never  given  as  functions  of  the  Richardson  number  and  surface  rough- 
ness, despite  extensive  micrometeorological  experience  that  this  is  so.   Never- 
theless, it  is  our  opinion  that  this  is  of  far  lessor  importance  than  the 
measurement  and  methodological  questions  to  be  discussed  next. 

Typical  measurement  programs  in  air  pollution  meteorology  are  grossly  in- 
adequate.  The  example  discussed  above  was  a  good  state-of-the-art  effort  which 
attempted  to  apply,  with  great  energy  and  professional  competence,  known  tech- 
niques and  standard  methods  to  the  problem.   Nevertheless,  it  is  clear  that  the 
information  and  results  obtained  are  not  adequate  for  the  purpose  of  making 
future  projections.   Thus,  a  general  and  serious  constraint  must  be  faced.   If 
we  wish  to  make  reliable  projections  of  the  air-quality  impact  of  highway 
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Figure  1-16.  Projections  to  the  year  1990  for  air  quality  with 
and  without  a  proposed  freeway. 
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development,  we  must  do  better  than  we  are  now. 

Specifically,  present  measurement  programs  have  the  following  faults. 
Air  concentration  of  pollutants,  diffusion,  and  standard  meteorological  par- 
ameters are  measured,  but  the  Richardson  number  is  not.   This  deficiency  en- 
sures lack  of  repeatability  in  the  results  and  makes  an  unacceptably  large 
scatter  in  the  data  inevitable.   The  Richardson  number  has  been  shown  to  be 
the  single  most  important  parameter  governing  all  micrometeorological  processes. 
It  should  be  a  standard  part  of  any  air-quality  measurement  program.   Secondly, 
measurements  are  not  adequately  extensive  either  in  time  or  space.   Two  and  a 
half  months  is  not  a  sufficiently  long  time  to  develop  reliable  meteorological 
or  air-quality  statistics.   Likewise  measurements  of  air  quality  at  various 
distances  from  a  highway  must  be  simultaneous  or,  again,  the  data  will  not  be 
interpretable . 

It  may  be  argued  that  the  measures  which  are  implied  in  the  criticism 
above,  i.e.,  more  sophisticated  and  extensive  measurements,  would  be  too  ex- 
pensive to  fund.   The  counter  argument  is  that  present  practices  are  not  pro- 
ducing the  results  we  need.   It  would  be  better  to  run  far  fewer  measurements 
on  a  higher  plane;   the  results  would  be  far  more  useful. 

Until  adequate  measurements  are  available,  the  question  of  model  valida- 
tion remains  academic.   However,  considerably  more  effort  could  and  should  be 
spent  in  analysis  of  the  properties  of  the  models  themselves.   Intensive  sen- 
sitivity studies  reveal  acceptable  and  unacceptable  properties  of  models.   In 
addition,  such  sensitivity  information  helps  put  the  instrumental  and  obser- 
vational effort  in  better  perspective.   For  instance,  it  doesn't  make  sense 
to  make  an  expensive  effort  to  measure  an  input  parameter  which  has  little 
effect  on  the  final  model  output. 
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An  Integrated  Microclimate/Diffusion  Approach 

In  the  previous  sections,  we  have  presented  a  survey  of  the  factors  which 
control  wind  transport,  diffusion  and  fallout  of  pollutants  emitted  by  motor 
vehicles.   In  this  section,  I  shall  describe  the  research  approach  of  the  at- 
mospheric environment  group  at  the  University  of  California,  Davis.   This 
effort  was  a  cooperative  effort  among  the  Departments  of  Agricultural  Engineer- 
ing, Geography,  Water  Science  and  Engineering,  Civil  Engineering,  the  Division 
of  Environmental  Studies  and  the  Atmospheric  Science  group.   The  work  I  shall 
be  outlining  here  is  largely  a  group  effort  and  my  colleagues,  Donald  Morgan, 
Wilson  Goddard,  Vinaya  Sharma ,  and  David  Rogers,  share  in  the  credit  for  our 
results. 

Our  overall  objective  was  to  study  the  fundamental  processes  which  give 
rise  to  the  diverse  microclimates  which  are  observed  in  urban  areas.   We  wished 
to  develop,  as  consequences  of  the  most  basic  theory  available,  means  for  cal- 
culating atmospheric  stability  and  diffusion  parameters  and  hence,  air  quality. 
One  final  objective  of  our  research  programs,  which  we  have  only  partially  a- 
chieved,  is  to  develop  models,  not  only  of  air  pollution  concentration,  but 
also  of  human  comfort  and  health. 

We  selected  the  city  of  Sacramento,  California  as  the  target  area  for  this 
study  because  of  topographic  simplicity,  relative  isolation  from  other  urban 
areas,  stable  summer  meteorological  conditions  and  proximity  to  the  Davis 
Campus.   Figure  1-17  shows  the  geographical  location  of  Sacramento  within  the 
Central  Valley  of  California,  and  Figure  1-18  shows  a  map  of  the  city  with  the 
grid  we  used  to  describe  the  properties  of  the  city.   Each  rectangle  measures 
6,000  to  8,000  feet  and  hence  the  area  is  48  million  square  feet  or  approximately 
1.7  square  miles. 

We  began  our  study  by  dividing  the  city  into  13  conventional  land-use 
categories.   Table  1-4  shows  the  13  categories  and  their  relative  abundance  in 
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Figure  1-17.  Wind  flow  in  the  Central  Valley  of 
California. 
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Figure  1-18.  The  City  of  Sc.cranento  showing  the  spatial  grid  used  in 
the  micaxxrliirate/diffusior  st-jdy. 
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TABLE  1-4. 

Category 

Abbreviation 

Light-Density  Residendial 

LDR 

Medium-Density  Residential 

MDR 

Heavy-Density  Residential 

HDR 

Schools 

SCH 

Shopping 

SHP 

Office  Buildings 

OFF 

Parks 

PRK 

Central  Business  District 

CBD 

Industrial 

IND 

Freeways 

FWY 

Open  Green 

PGR 

Seasonal  Green 

SGR 

Water 

WTR 

%   of   City 
43.6 

1.3 

0.6 

3.2 

8.7 

0.6 

8.2 

0.3 
15.0 

2.7 

1.5 
11.8 

2.5 
100.0 
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our  target  area  as  a  whole.   As  can  be  seen,  the  categories  of  Light-Density 
Residential,  Industrial,  and  Seasonal  Green  (agricultural  land)  occupy  over 
70%  of  the  study  area.   From  either  the  known  properties  of  each  land-use  cate- 
gory or  field  measurements  of  our  own,  we  were  able  to  attach  numbers  to  each 
category  describing  its  fundamental  physical  properties.   Table  1-5  shows  the 
result  of  this  effort.   The  aerodynamic  roughness  length,  z  ,  was  obtained  by 
means  of  a  technique  due  to  Lettau  (1969)  in  which  geometric  measurements  of 
the  average  properties  of  the  roughness  elements  (buildings  and  trees)  can  be 
used  to  infer  this  parameter.   The  "wet  fraction"  is,  in  fact,  the  percentage 
area  occupied  by  green,  transpiring  plants  in  a  given  land-use  category. 

All  micrometeorological  experience  (Sutton, 1953 ;  Sellers, 1965  ;  Munn,1966) 
indicates  that  these  are  the  fundamental  physical  and  biological  characteristics 
of  the  surface  of  the  earth  which  determine  the  "response"  of  a  given  area  to 
large-scale  meteorological  processes,  such  as  radiation  load,  wind  speed,  or 
temperature  and  humidity  well  away  from  the  surface.   The  "response"  is,  in 
fact,  the  microclimate  which  determines,  in  turn,  atmospheric  stability,  dif- 
fusion parameters,  and  human  comfort.   Given  detailed  information  about  the 
nature  of  the  surface  of  the  earth,  it  is  possible  to  calculate  microclimatic 
parameters.   The  heat-budget  equation  for  the  surface  of  the  earth, 

R=  H  +  LE  +  G  ,(20) 

which  is  an  expression  of  the  first  law  of  thermodynamics,  is  the  single  most 
powerful  principle  which  determines  the  microclimate.   In  equation  (20),  R^  is 
the  net  flux  of  radiation  absorbed  at  the  surface  of  the  earth  (positive  dur- 
ing the  day,  negative  at  night),  H  is  the  flux  of  heat  carried  upward  from  the 
warm  surface  of  the  earth  by  turbulent  convection,  E  is  the  evaporation  rate 
(including  transpiration) ,  and  L  is  the  latent  heat  of  evaporation,  so  that 
the  product,  LE ,  is  the  flux  of  energy  required  to  allow  evaporation  rate,  E, 
and  G  is  the  flux  of  heat  into  the  ground.   The  net  radiation  term  can  be 
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TABLE  1-5.   CITY  INTERFACE  PARAMETERS  BY  LAND-USE  CATEGORY 


Albedo   Infrared  Heat  Con-  Heat  Ca-   Wet 

Emissivity  ductivity  pacity   Fraction 

[cal  cm  [calcm- 
sec-1^"1]   "C"1] 


Roughness, z 
[cm] 


Light-density 
Residential 


.20 


88 


0023 


.45 


.43 


108 


Medium-density 
Residential 

.23 

.86 

.0024 

.43 

.56 

532 

Heavy-dens  ity 
Residential 

.25 

.88 

.0022 

.40 

.43 

370 

Schools 

.15 

.90 

.0023 

.50 

.35 

41 

Shopping 

.20 

.86 

.0022 

.42 

.04 

36 

Office 

.22 

.85 

.0025 

.45 

.13 

175 

Park 

.16 

.89 

.0030 

.58 

.93 

127 

Central  Business 
District 

.26 

.86 

.0024 

.42 

.07 

321 

Industrial 

.26 

.86 

.0024 

.42 

.01 

13 

Freeway 

.30 

.80 

.0030 

.60 

.00 

4 

Open  Green 

.35 

.90 

.0006 

.30 

.05 

2 

Seasonal  Green 

.15 

.90 

.0030 

.60 

1.00 

2 

Water 

.09 

.93 

.0014 

1.00 

1.00 

.001 
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thought  of  in  simplified  terms  as  the  energy  input  rate  while  the  terms  on  the 
right  determine  how  the  energy  is  used.   For  instance,  over  water  or  surfaces 
completely  covered  with  freely  transpiring  plants,  the  majority  of  the  incoming 
solar  energy  goes  into  evaporating  water  (LE)  and,  hence,  is  not  available  to 
heat  the  air  or  underlying  surface.   On  the  other  hand,  for  surfaces  which  are 
dry  and  do  not  evaporate,  most  of  the  sun's  energy  goes  into  heating  the  air 
(H)  or  surface  material  (G) .   This  simple  mechanism  explains  most  of  the  micro- 
climatic variability  found  in  mesoscale  systems  such  as  the  urban  environment. 
A  glance  at  our  "wet  fraction"  data  shows  that  very  large  contrasts  exist  be- 
tween the  various  land-use  categories  in  the  amount  of  evaporation  which  must 
be  occurring.   This  also  is  largely  the  basis  for  contrasts  in  stability,  dif- 
fusion, air  quality  and  human  comfort. 

Solution  of  the  energy  budget  equation  is  feasible  using  micrometeorological 
formulae  for  each  of  the  various  terms.   The  technique  has  been  successfully 
used  by  Halstead  (1957) ,  Myrup  (1969) ,  Goddard  (1971) ,  and  Outcalt  (1972) .   The 
output  of  these  models  includes  air  and  surface  temperature  as  well  as  the 
energy-budget  components  and  the  friction  velocity,  UA ,  discussed  earlier.   It 
should  be  remembered  that  the  only  two  important  variables  which  enter  into  the 
definition  of  the  Monin-Obukhov  length,  L,  are  H  and  u\.   Hence,  we  are  in  a 
position  to  calculate  the  microclimate  and  atmospheric  stability  from  interface 
characteristics . 

We  have  performed  this  calculation,  using  the  microclimate  simulation  model, 
CLISIM,  developed  by  Goddard  (1971).   CLISIM  requires,  as  input,  the  interface 
parameters  as  listed  above,  meteorological  parameters  (temperature,  wind  speed, 
and  humidity)  well  away  from  the  surface,  and  latitude  and  date  (to  determine 
the  solar  radiation) . 

As  a  preliminary  feasibility  study,  we  have  performed  the  CLISIM  calcula- 
tions for  each  of  our  13  land-use  categories  for  August  24,  1971.   Table  1-6 
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TABLE  1-6.   FRICTION  VELOCITY  AND  MONIN-OBUKHOV  LENGTH 
BY  LAND-USE  CATEGORIES  FOR  AUGUST  ]4,  1971 


U^(cm 

s"1) 

0500 

1400 

LDR 

29.7 

53.1 

MDR 

38.7 

60.5 

HDR 

37.5 

59.3 

SCH 

15.7 

47.7 

SHP 

— 

55.1 

OFF 

41.4 

61.8 

PRK 

14.4 

44.6 

CBD 

40.0 

69.8 

IND 

— 

49.0 

FWY 

5.8 

44.9 

OGR 

— 

42.5 

SGR 

5.9 

31.9 

WTR 

9.7 

26.5 

L(met< 

ers) 

0500 

1400 

56.29 

-  55.01 

50.60 

-243.40 

58.49 

-127.00 

12.20 

-  40.00 

— 

-  29.70 

84.53 

-  52.71 

9.04 

-139.49 

94.62 

-  62.26 

-  24.56 

3.06 

-  19.64 

— 

-  17.69 

3.26 

-  27.17 

•  20.94 

-  86.50 
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shows  the  result  of  this  operation  for  the  friction  velocity,  U^ ,  and  the 
Monin-Obukhov  length,  L,  at  0500  and  1400  PST.   Pronounced  variation  among 
the  various  land-use  categories  is  apparent. 

The  real  city  is  not  made  up  of  neatly  arranged  parcels  of  pure  land-use 
categories.   Instead  every  one  of  our  138  1.7-square  mile  zones  contained  sev- 
eral categories.   We  have  measured  the  percent  of  each  land-use  category  in 
each  rectangle.   With  this  information,  we  can  form  a  weighted  average  for  the 
pure  land-use  parameters  and  produce  a  number  which  should  be  representative 
of  the  entire  1.7-square  mile  area.   Figure  1-19  shows  the  result  of  such  an 
analysis  for  the  roughness  parameter,  z  .   The  city  appears  as  a  rough  oasis 
in  the  midst  of  relatively  smooth  agricultural  lands.   Z  values  are  partic- 
ularly big  in  the  vicinity  of  the  central  business  district  due  to  the  presence 
of  many  large  buildings.   We  have  made  similar  analyses  of  all  the  interface 
parameters  listed  above  and  for  the  energy-budget  components.   Figure  1-20 
shows  the  result  for  R^j  net  radiation,  and  Figure  1-21  for  the  evaporation 
term,  LE,  for  the  1400  PST  time  period.   At  this  time,  the  city  is  absorbing 
less  radiation  than  the  surrounding  agricultural  lands  due  to  higher  reflectiv- 
ity and  greater  emission  of  infrared  energy.   In  the  evaporation  display,  the 
city  appears  as  a  dry  island  surrounded  by  large  evaporation  from  the  agri- 
cultural areas. 

We  have  now  reached  the  point  where  we  have  developed  the  capability  of 
calculating  the  atmospheric  stability  from  measurable  surface  properties,  in- 
cluding its  roughness.   The  next  step  would  be  to  specify  the  air  pollution 
emission  process  and  finally  the  transport  and  dispersion.   This  we  have  done, 
and  we  briefly  describe  our  approach  in  the  remainder  of  this  section.   Figure 
1-22  shows  the  emission  of  carbon  monoxide  from  motor  vehicles  in  Kilograms 
per  day  per  1.7  square  miles.   These  figures  were  estimated  as  follows.   For 
each  rectangle,  the  total  miles  of  residential,  arterial,  and  highways  were 
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Figure  1-19.  Aerodynamic  roughness 
parameter,  Z„. 
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Figure  1-20.  Net  radiation  flux  for  1400 
PST,  August  24,  1971  according  to  CLISIM. 
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Figure  1-22.  Total  emission  of  carton 
monoxide  per  24  hours  per  1.7  square 
miles. 
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measured  from  maps.   The  average  traffic  volume  for  each  of  the  three  types  of 
streets  was  obtained  from  the  California  Division  of  Highways  and  car-mile 
factors  computed.   Application  of  vehicle  emission  factors  (U.  S.  Env.  Pro- 
tection Agency,  1972)  at  appropriate  average  vehicle  speeds,  then  yields  the 
total  amount  of  each  kind  of  pollution  emitted  by  the  three  kinds  of  streets, 
and,  after  summation,  the  total  emitted  by  the  rectangle  in  question. 

The  pattern  of  carbon  monoxide  emission  from  motor  vehicles  varies  im- 
mensely from  area  to  area.   The  traffic  concentration  observed  in  the  central 
business  district  results  in  the  largest  value  found,  21,258  Kilograms  per  day. 
At  the  edge,  in  some  agricultural  areas,  we  show  zero  values  although  this  is 
undoubtedly  an  underestimate  if  only  because  we  did  not  take  agricultural  ve- 
hicles into  consideration. 

We  have  developed  a  preliminary  air  pollution  transport  and  diffusion 
model  which  is  designed  to  take,  as  input,  the  kind  of  information  which  we 
have  developed  for  the  state  area.   The  model,  which  was  the  subject  of  a  por- 
tion of  Dr.  Sharma's  dissertation,  is  related  to  the  work  of  Monin  (1959)  and 
can  be  considered  to  be  a  member  of  the  similarity  family  of  diffusion  models. 
It  presently  exists  in  a  form  which  neglects  the  effects  of  elevated  inversion 
layers  and  has  been  formulated  only  for  an  infinite  line  source  configuration. 
These  restrictions,  however,  are  not  inherent  in  the  approach,  and  extension 
of  the  model  to  include  more  realistic  features  is  planned. 

The  central  physical  phenomenon  of  interest  is  the  behavior  of  the  plume 
of  pollutant  which  is  emitted  by  a  highway  line  source.   The  vertical  pro- 
pagation by  turbulent  diffusion  and  horizontal  transport  of  the  plume  determines 
the  resultant  air  concentration  of  pollutant  observed  at  some  downwind  location. 
Figure  1-23  shows  the  model  geometry.   Monin  (1959)  postulated  that  "turbulent 
diffusion  in  a  horizontally  homogeneous,  stationarysurface  layer  of  air  obeys 
the  similarity  theory  in  which  L  and  U^  are  the  only  scales  of  length  and 
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Plume  boundary 


Mean  wind  u(z) 
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^Ground  (complete  reflector) 


Infinite  line  source 
of  strength  =  Q 


Figure  1-23.  Problem  description. 
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velocity".   It  follows,  then,  from  dimensional  considerations  that  the  upward 
velocity  of  propagation  of  marked  or  polluted  air,  due  to  diffusion,  W^ ,  is 
given  by 

W^  =  A  U^  4>  (z/L)  ,(21) 

where  U^  and  L  are  defined  as  before,  A  is  a  constant  with  the  value  of  W^/U^ 
for  neutral  stratification  (z/L  =  0).   The  equation  which  governs  the  height, 

Z,  of  the  upper  boundary  of  plume,  is 

W 

^  =  —  (22) 

dx   U  lK      } 

where  x  is  downwind  distance.   Monin  assumed  the  "log-linear"  relationship  for 
U  (a  stability-corrected  form  of  the  logarithmic  wind  law)  and  obtained 

dz  .    k  A  4>  (z/L) 


dx 


(ln(z/z  )  +  6  (—2-) 

O  L 


,(23) 


where  $  is  another  constant.   From  this  formulation,  Monin  further  postulated 
that  "The  shape  of  the  boundaries  of  the  plume  (in  particular,  their  inclina- 
tion to  the  horizontal)  does  not  depend  upon  the  wind  velocity  but  does  depend 
upon  the  stratification  of  the  atmosphere".   Monin  went  on  to  integrate  the 
above  equation  using  empirical  information  for  the  ^-function  and  claimed  good 
correspondence  between  calculations  and  smoke-plume  observations.   We  have  fol- 
lowed the  general  principles  proposed  by  Monin  but  have  chosen  to  use  different 
techniques  to  specify  the  plume  shape  and  velocity  profiles. 

The  model  is  constrained  to  obey  the  principle  of  conservation  of  mass  as 

stated  by 

z 

Q  =  /  U(z)  C(z)  dz  ,(24) 

o 

where  Q  is  the  source  strength.   For  mathematical  convenience,  and  also  in 
hopes  of  validity  at  greater  heights,  we  represent  the  wind  by  a  power  law  form 

U  -  IL  (-  )P  ,(25) 

where  \J     is  the  wind  speed  at  height  z  .   The  exponent  must  be  a  function  of 
surface  roughness  and  stability.   Davenport  (1965)  gives  empirical  values  for  p 
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which  range  from  0.1  for  open  sea,  to  0.16  for  flat  open  country  to  0.28  for 
forest  to  0.4  for  urban  areas.   Reliable  information  on  the  stability  dependence 
of  p  is  not  available.   In  the  absence  of  such  information,  we  were  forced  to 
use  the  log-linear  formula  to  generate  p-values  as  functions  of  both  stability 
and  roughness.   This  procedure  is  only  an  interim  measure  until  better  data 
are  available.   The  values  of  p,  so  obtained,  do  have  the  virtue  of  allowing 
convenient  mathematical  manipulations  in  the  model  and  yield  reasonable  ex- 
trapolations to  heights  where  the  log  +  linear  form  alone  fails.   Figure  1-24 
shows  p-values  for  stable  stratification  as  a  function  of  z   and  L,  and  Figure 
1-25  shows  the  same  display  for  the  unstable  case.   Most  surprisingly,  these 
values  correspond  quite  well  to  Davenport's  values. 

In  our  approach  we  obtain  the  height  of  the  top  of  the  plume  by  means  of 
an  analytical  solution,  due  to  Roberts  (1923),  of  the  plume-diffusion  equations 
with  wind  speed  and  diffusivity  varying  according  to  a  constant  power  of  height 
(power  law) .   After  considerable  manipulation,  a  result  is  obtained  for  the 
height,  z,  of  the  plume 

z  =  a  x  ,  (26) 

where  b  =  l/(2p  +  1)  and  a  is  a  known  function  of  stability  and  roughness. 
Figure  1-26  shows  an  example  of  the  growth  of  a  polluted  plume  over  an  area  of 
large  roughness  (z  =  100  cm)  for  various  stabilities. 

The  air  concentration  observed  at  any  point  within  the  plume  depends  on 
the  form  of  the  vertical  distribution  of  pollutant.   The  available  data  are 
not  particularly  helpful  in  this  respect.   Figure  1-27  shows  measurements  of 
smoke  density  from  Monin  (1959) ,and  Figure  1-28  data  from  the  Fort  Wayne  ex- 
periment conducted  by  the  Traveler's  Research  Center  (Hilst  and  Bowne,  1966). 
In  both  cases,  and  in  all  others  we  have  inspected,  the  data  could  be  repres- 
ented by  a  linear  decrease  with  height  as  well  as  any  other  function.   We  have 
chosen,  therefore,  to  assume  a  linear  decrease  of  concentration  with  height 
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Figure  1-24.  The  exponent  in  the  power  law,  p,  vs.  the  stratification 
parameter,  L,  for  different  roughnesses,  Ze,  for  stable  stratification. 
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Figure  1-26.  Height  of  the  plune,  H,  vs.  downwind  distance,  X 
for  surface  roughness,  Z0  =  1.0  m. 
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Figure  1-27.  The  vertical  distribution  of  smoke  density 
in  the  smoke  plume,  from  Monin  (1959) . 
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Figure  1-28.     Vertical  dosage  distribution  at  WWEyiV 
tcwer,  first  sanpling  line,  Ft.  Wayne  "check"  trial 
THC  65-16-G1.     Fran  Hilst  and  Bcwne  (1966). 
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within  the  plume.   This  will  probably  lead  to  incorrect  concentrations  near 
the  upper  edge  of  the  plume  since  the  concentration  must  assume  a  more  asym- 
ptotic shape  there.   However,  the  linear  representation  of  concentration 
throughout  the  majority  of  the  plume  should  be  as  good  as  the  observational 
data.   Integration  of  the  conservation  of  mass  equation  using  the  assumed  wind 
profile,  plume  height,  and  concentration  profile  formulae  yields  a  relationship 
for  the  surface  concentration,  C  , 

c  -  _2_  (%v    (p+2)  (p+1) (27) 

o  "  UXZ  ^z  '    (p+2)  +  a^p+1)  tK    '} 


wher 


e  U1  is  the  wind  speed  at  height  z,  and  a.,  is  constant, 


We  have  thus  obtained   this,  a  simple  diffusion  and  transport  model  which 
requires ,  as  input,  surface  roughness,  stability,  wind  speed  at  a  specified 
height  and  the  source  strength.   We  feel  that  this  model  is  probably  as  good 
as  the  data  and  has  the  virtue  of  being  explicitly  dependent  on  the  most  im- 
portant physical  parameters  which  are  known  to  govern  the  diffusion  process 
in  the  atmosphere. 

In  order  to  analyze  the  sensitivity  of  the  model  to  the  various  para- 
meters which  determine  the  surface  concentration,  we  define  the  sensitivity 

of  the  model  to  parameter  s,  Ss,  as 

C 
100   3  o 
S  =  ~-  - —  •  (28) 

s    C    8s  v  ' 

o 

The  algebraic  simplicity  of  this  model  allows  immediate  analysis  of  the 
sensitivity  to  wind  speed,  U,  source  strength,  Q,  and  plume  height,  H.   These 
sensitivities  are  given  below. 

u  "    U 


S  =   100 


sH  =  -100  (P+1) 


H  2p  +3  (29) 
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The  model  is  most  sensitive  to  parameter  change  or  error  at  low  wind  speed, 
weak  source  strength,  and  low  plume  heights. 

Sensitivity  of  the  model  to  surface  roughness  and  stability  was  estimated 
by  performing  calculations  for  specific  conditions.   For  instance,  the  rough- 
ness calculation  was  done  at  a  distance  of  3  Km  downwind  for  neutral  strati- 
fication and  for  roughnesses  of  10  and  200  cm.   Under  these  conditions,  the 
prediction  was  that  the  smooth  surface  concentration  would  be  2.4  times  that  of 
the  rough  surface  case.   A  similar  strategy  was  followed  to  obtain  the  sen- 
sitivity to  stability.   For  this  calculation,  the  roughness  parameter  was  set 
at  200  cm,  a  wind  speed  of  1  ms   at  5  meters,  and  an  arbitrary  source  strength 
of  1000,  and  we  obtained  the  surface  concentrations  given  in  Table  1-7.   The 
concentration  calculated  at  L  =  100  (stable)  is  5.5  times  that  at  L  =  -100  (un- 
stable) . 

Direct  verification  of  this  or  any  other  model  is  not  straight  forward 
because  of  the  lack  of  definitive  experimental  evidence.   For  instance,  the 
Fort  Wayne  study,  reported  by  Hilst  and  Bowne  (1966),  is  by  far  the  best  avail- 
able set  of  data.   In  this  experiment  a  long  line  of  uniform  fluorescent  par- 
ticles was  released  upwind  of  the  city  by  airplane  and  resultant  air  concen- 
trations were  measured  in  the  city  at  five  points  as  well  as  stability  and 
turbulence  parameters.   Unfortunately,  the  source  strength  was  not  well  known 
in  this  experiment.   The  manner  in  which  the  particles  are  generated  and  re- 
leased in  the  wake  of  the  airplane  results  in  the  immediate  deposition  of  a 
large  fraction  of  the  particles.   For  this  reason  Hilst  and  Bowne  recommend, 
for  the  purposes  of  verification  calculations,  that  the  source  strength  be 
adjusted  to  yield  observed  values  at  the  first  sampling  point  downwind.   Using 
observed  wind  and  stability  (L  =  -200  meters)  and  a  roughness  of  200  cm,  we 
obtained  the  results  given  in  Table  1-8  in  terms  of  dossage,  D  .   With  the  ad- 
justment at  1600  meters,  the  correspondence  between  observed  and  calculated 
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TABLE  1-7. 

Monin  -  Obukhov  Calculated  surface 

Length, L(ra)  concentration^ 


100  0.96 

500  2.44 

-10,000  3.12 

+10,000  3.38 

+   500  3.86 

+   100  5.32 
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TABLE  1-8. 
OBSERVED  AND  CALCULATED  SURFACE  DOSSAGE  FOR  FORT  WAYNE  TEST  NO. 65-06  Gl 


Ratio 


X 

Observed 

Calculated 

(Km) 

D   (particle- 
o 

-min/m  ) 

1,600 

96.2 

96.2 

4,800 

42.0 

48.6 

8,000 

35.3 

32.4 

11,200 

31.5 

24.8 

14,400 

24.5 

21.2 

1.0 

0.87 
1.09 
1.27 
1.15 


TABLE  1-9. 
OBSERVED  AND  CALCULATED  DOSSAGE  FOR  FORT  WAYNE  TEST  NO. 65-13  Gl 

Ratio 


1.0 

0.84 

0.87 

0.74 

0.78 


X 

Observed 

Calculated 

(Km) 

D   (particle- 
o 

-min/m  ) 

3,200 

54.5 

54.5 

6,400 

36.2 

43.2 

9,600 

31.2 

36.0 

12,800 

23.0 

31.0 

16,000 

21.6 

27.0 
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dossages  is  good.   For  these  conditions,  the  model  generated  an  exponent  in 
the  wind  speed  power  law  of  0.30;   the  observed  value  from  tower  measurements 
in  the  city  was  0.32. 

A  companion  calculation  was  made  for  a  slightly  stable  case  in  which  L 
was  measured  at  510  m.   Again,  adjusting  the  source  strength  to  force  exact 
correspondence  at  the  first  downwind  measurement  point  and  using  the  observed 
wind  speed,  we  obtained  the  results  shown  in  Table  1-9.   In  this  case  the  model 
appears  to  over-predict  air  concentrations  by  about  20%  at  considerable  downwind 
distances. 

This  brief  verification  study  is  meant  to  be  only  illustrative.   The 
available  experimental  data  do  not  permit  unambiguous  conclusions  to  be  drawn 
about  the  limits  of  accuracy  of  any  model.   The  data  do  suggest  that  our  model 
may  be  a  useful  addition  to  the  tools  we  have  at  our  disposal  to  predict  and, 
ultimately,  control  air  pollution.   Our  goal  was  to  formulate  a  simple,  easy 
to  use,  model  which  required  surface  roughness  and  stability  as  primary  inputs 
and  which  yielded  results  comparable  to  the  quality  of  the  data  available. 

In  particular,  we  wished  to  develop  means  of  handling  the  dispersion 
processes  which  were  closely  related  to  the  primary  influence  on  turbulence 
near  the  surface,  i.e.,  roughness  and  stability.   It  appears  that  we  may  have 
come  close  to  achieving  our  goals. 

One  of  our  goals  is  to  provide  information,  conceptual  and  analytical 
tools,  which  will  be  useful  in  the  decision-making  and  planning  portion  of 
society.   Figure  1-1  suggests  and  certainly  our  analysis  strongly  supports  the 
idea  that  land-use  is  a  particularly  critical  portion  of  the  processes  that  air 
pollution  is  a  part  of.   Ideally,  then,  we  would  hope  to  relate  pollution  emis- 
sion and  air  quality  to  land-use  categories.   We  are  in  the  process  of  doing 
this.   Table  1-10  presents  a  preliminary  analysis  of  the  emission  of  carbon 
monoxide  from  motor  vehicles  by  land-use  type.   Emissions  are  presented  as  an 
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037 

.177 

0 

.203 

026 

1.511 

006 

.464 

019 

.434 

0 

.019 

010 

.405 

002 

.413 

024 

.476 

095 

.298 

0 

1.537 

165 

.635 

339 

.792 

0 

.019 

0 

.415 

036 

.451 

2 
emission  flux  of  grams  of  CO  per  ft  .   The  data  were  selected  from  spot  read- 
ings of  representative  portions  of  our  spatial  array. 

TABLE  1-10.   PRELIMINARY  ANALYSIS  OF  CARBON  MONOXIDE  EMISSION  FROM 

MOTOR  VEHICLES  BY  LAND-USE  TYPE  FOR  SACRAMENTO,  CALIFORNIA 

2 
in  grams  CO/ (feet) 

Category Residential Arterial Highway Total 

LDR 

MDR 

HDR 

SHP 

IND 

OFF 

CBD 

SGR 

We  expect  the  general  features  of  this  table  to  remain  the  same  when  the 
analysis  is  completed.   In  particular,  the  fact  that  the  emission  from  arter- 
ial streets  is  larger  than  highway  has  been  a  persistent  feature  of  our  pre- 
liminary results.   Largest  emission  fluxes  are  found  in  the  Heavy-Density 
Residential,  Industrial,  and  Shopping  land-use  categories. 

It  may  be  that  emission  factors  like  these  may  be  typical  of  similar  land- 
use  categories  throughout  California.   We  plan  to  pursue  this  concept  in  our 
future  research.   We  hope  that  future  planning  decisions  regarding  land-use 
will  take  into  account  probable  emission  of  pollutants  and  resultant  impact  on 
air  quality. 

To  close  this  section,  we  present  a  second  preliminary  result  of  our  systems 
approach  to  air  quality.   The  final  result  of  the  current  phase  of  our  research 
program  has  been  a  data  base  consisting  of  all  the  information  necessary  to 
make  an  air  pollution  transport  and  dispersion  calculation  with  the  Sharma  model, 
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presented  in  the  previous  section.   In  order  to  provide  an  example  of  the  use 
of  air  quality  modeling,  we  have  performed  a  preliminary  calculation  using 
the  Sacramento  data  base  and  CLISIM  estimates  of  stability  as  inputs  to  the 
Sharma  model.   To  do  this,  we  assumed  an  overall  value  of  66  cm  for  the  rough- 
ness parameter  and  the  following  values  of  CO  emission  by  18,000-foot  strips 
(the  city  imagined  infinite  for  the  purpose  of  this  calculation  in  the  cross- 
wind  direction) . 

TABLE  1-11.   CO  EMISSION  AS  A  FUNCTION  OF  DOWNWIND  DISTANCE 

-2    -1 
in  gm  CO  m  sec 

X(f t)  Area  Emission  Rate 

18,000  9.4  x  10~6 

36,000  16.8  x  10"6 

54,000  11.7  x  10"6 

72,000  9.6  x  10"6 

90,000  6.2  x  10"6 

108,000  2.5  x  10~6 

These  emission  rates  were  chosen  to  be  representative  of  the  pattern  in  Sacra- 
mento during  the  normal  summertime  wind  flow  from  the  Southwest.   We  then  per- 
formed an  air  quality  calculation  for  the  times  of  0500  and  1400  PST  using 
average  L-values  generated  by  CLISIM  for  the  August  24,  1971  case.   These  were 
15  meters  for  the  stable  case  and  -40  meters  for  the  unstable  afternoon  case. 

The  results  of  this  exercise,  which  is  meant  to  be  only  illustrative,  are  shown 

-3   3 
in  Figure  1-29  and  Figure  1-30.   The  CO  concentration  is  given  in  gm  m  xlO  . 

The  values  on  the  vertical  ordinate  can  be  taken  as  ppm  for  sea  level  conditions 

2 
(ignoring  the  factor  of  10  ) .   The  unstable  case  shows  a  pattern  of  CO  concen- 
tration which  follows  the  assumed  emission  rate  very  closely.   Concentrations 
are  low.   Both  these  results  are  the  consequence  of  the  efficient  dispersion 
characteristic  of  unstable  stratification.   In  contrast,  the  values  of  concen- 
tration calculated  for  the  stable  0500  PST  case  are  very  high  and  the  shape  of 
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Figure  1-30.  Confuted  Carbon  Monoxide  concentration. 
Stable  case. 
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Figure  1-29.  Canpated  Carbon  Monoxide  concentration. 
Unstable  case. 
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the  concentration  curve  is  different.   The  effect  of  the  very  stable  atmosphere 
is  a  downwind  buildup  of  air  concentration  to  levels  about  35  times  that  cal- 
culated for  unstable  conditions.   It  should  be  noted  that  we  used  the  same 
traffic  density  for  both  of  these  order-of-magnitude  calculations. 
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Summary  and  Conclusions 

In  this  report,  we  have  attempted  to  tie  together  the  physical  processes 
which  generate  the  microclimate,  determine  transport  and  diffusion  processes 
on  the  micro-  and  meso-scale  and  how  these  are  related  to  land-use  patterns 
and,  finally,  air  quality.   There  are  four  general  areas  in  which  conclusions 
from  this  review  are  appropriate. 

(1)  The  microclimate  emerges  from  this  discussion  as  a  concept  of  con- 
siderable importance  to  the  problem  of  measuring  and  predicting  air  quality. 
Large  horizontal  variations  of  important  physical  processes,  such  as  those  re- 
presented by  the  flux  of  net  radiation,  exist  within  cities.   These  "neighbor- 
hood contrasts"  are  strongly  coupled  to  land-use  patterns  and  hence  directly 
related  to  the  planning  processes  in  urban  and  regional  governments.   In  par- 
ticular, the  evaporation  rate,  as  indicated  by  relative  abundance  of  green, 
freely  transpiring  plant,  is  the  single  most  important  parameter  which  de- 
termines the  microclimatic  response  of  a  specific  locality  to  a  given  radiation 
load  and  other  large-scale  meteorologic  factors.   Since  air  pollution  emission 
from  motor  vehicles  is  also  closely  related  to  land-use  patterns,  the  "air  pol- 
lution system"  as  a  whole  is  thus  strongly  coupled  to  land-use. 

(2)  The  processes  through  which  pollution  diffuses  from  highways  are  also 
largely  determined  by  the  microclimate.   In  particular,  surface  roughness  and 
atmospheric  stability,  as  measured  by  the  Richardson  number  or  its  equivalent, 
the  Monin-Obukhov  length,  are  the  most  important  parameters  controlling  micro- 
scale  transport  and  diffusion.   In  view  of  its  central  importance  to  micro- 
meteorological  and  diffusion  processes,  an  important  conclusion  to  be  drawn 
from  this  review  is  that  measurement  of  the  Richardson  number  should  be  a 
standard  part  of  highway  air  quality  programs.   Extensive  micrometeorological 
experience  indicates  that  this  is  necessary  to  allow  measurements  made  at  dif- 
ferent locations  or  at  different  times  to  be  compared. 
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(3)  Consideration  of  the  multitude  of  air  pollution  models  which  are 
available  to  estimate  air  quality  leads  to  the  conclusion  that  the  models 
generate  far  more  information  at  greater  precision  than  is  available  from  field 
measurements.   In  fact,  we  feel  the  need  for  better  and  more  comprehensive 
measurements  is  the  most  important  conclusion  to  be  drawn  from  this  review. 

The  measurement  problem  is  severe  because  large  local  variability  in  the  under- 
lying microclimate  leads  to  local  contrasts  in  all  processes  governed  by  the 
microclimate.   Hence,  any  measurement  program  must  be  carefully  designed  to 
ensure  that  the  measurements  are,  in  fact,  representative  of  the  area  in 
question.   In  particular,  pollution  and  meteorological  parameters  must  be  mea- 
sured simultaneously  within  the  study  region  over  an  adequate  averaging  time. 

(4)  We  have  also  briefly  reviewed  the  state  of  the  art  of  air  quality 
modeling  in  this  report.   The  prevailing  standard  in  the  field,  the  Gaussian 
family  of  models,  are  deficient  in  the  sense  that  the  diffusion  parameters 
which  are  used  in  these  models  are  not  derived  directly  from  the  physical  pro- 
cesses which  are  known  to  dominate  turbulent  diffusion  on  the  microscale,  i.e., 
the  microclimate  -  surface  roughness  -  Richardson  number  parameter  complex. 
More  generally,  it  appears  that  far  more  attention  needs  to  be  paid  to  the  pro- 
perties of  the  models  themselves.   Comprehensive  sensitivity  analyses  should 

be  made  in  conjunction  with  model  verification  programs.   Such  information 
will  allow  better  information  as  to  the  predictability  of  air  quality  and  of 
the  requirements  of  measurement  and  verification  programs. 
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VOLUME  III 
CHAPTER  2:  AIR  POLLUTION  CLIMATOLOGY 
FOR  SACRAMENTO,  CALIFORNIA 

Objectives 

As  a  portion  of  our  overall  program  of  research  into  the  factors  which  influence  the  at- 
mospheric environment  over  urban  areas,  we  have  studied  the  available  air  quality,  and  as- 
sociated meteorological  statistics,  for  Sacramento,  California,  for  the  period  of  1961-1969. 
We  had  three  general  objectives  for  this  work. 

(1)  During  the  period  1960-1970  the  population  of  Sacramento  County  increased  by  22% 
and  energy  consumption,  in  common  with  the  state  as  a  whole,  approximately  doubled.   In  very 
general  terms  it  might  be  expected  that  phenomena  associated  with  waste  and  inefficiencies, 
such  as  air  pollution,  would  show  increases  somewhere  between  the  rates  for  population  and 
energy  consumption,  other  factors  being  the  same.   This  last  condition  is,  of  course,  never 
met  in  practice.   The  major  meteorological  factors  which  determine  the  air  concentration  re- 
sulting from  a  given  rate  of  emission  of  pollution,  such  as  wind  speed  and  mixing  depth,  are 
not  constant  in  time.   In  addition,  the  technological  effort  to  control  automobile  emissions 
of  pollution  began  midway  through  the  1960-70  decade.   An  important  practical  objective  in 
examining  air  quality  data  for  this  period  would  be  to  determine  whether  a  measurable  im- 
provement can  be  demonstrated  from  these  data.   Therefore,  one  goal  for  this  study  was  to 
determine  whether  or  not  a  statistically  significant  trend  can  be  found  in  the  available  air 
quality  data.   We  recognize,  of  course,  that  an  air  pollution  "climatology"  may  not  even  be 
a  definable  concept  due  to  rapid  changes  with  time  and  lack  of  data. 

(2)  One  of  the  uses  to  which  historical  air  pollution  and  meteorological  data  may  be 
put  is  to  provide  baseline  information  for  model  development.   At  a  minimum,  the  historical 
data  sets  the  order  of  magnitude  of  general  pollution  levels,  as  well  as  maxima  and  minima, 
which  any  useful  predictive  model  must  be  able  to  generate  from  historical  emission  rates. 
In  addition,  it  would  be  hoped  that  sufficient  historical  data  are  available  that  detailed 
verification  tests  could  be  made  for  various  air  pollution  models.   Therefore,  the  second 
ovjective  of  our  study  of  historical  air  pollution  data  was  to  compile  sufficient  information 
to  test  predictive  models  of  the  processes  by  which  pollution  is  transported  and  diffused  by 
the  wind. 

(3)  As  a  third  objective,  we  wished  to  subject  the  available  data  to  standard  statistical 
procedures  to  determine  the  existence  of  strong  statistical  relationships  between  the  various 
pollution  and  meteorological  parameters.   Such  relationships,  if  found  at  a  sufficiently  high 
level  of  significance,  might  be  useful  for  various  pragmatic  applications  as  well  as  a  guide 
to  the  model-generating  process. 

Pollution  Data 

Data  for  contaminants  was  available  for  Sacramento  from  the  Air  Resources  Board  from 
October  1,  1957  to  December  31,  1969  at  the  time  we  initiated  our  study.   The  early  data  in- 
cludes only  oxidant  with  other  pollutants  appearing  gradually  as  appropriate  instrumentation 
became  available.   Pollution  data  for  the  entire  period  mentioned  above  were  punched  up  for 
computer  processing.   However,  since  meteorological  data  were  conveniently  available  only  for 
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January  1,  1961  to  December  31,  1969,  all  statistics  regarding  relationships  between  pollutant 
concentrations  and  meteorological  parameters  were  confined  to  this  shorter  period.   Data  were 
collected  for  every  day  in  the  target  period  for  the  following  parameters:  maximum  daily  con- 
centrations of  particulates,  total  hydrocarbons,  oxides  of  nitrogen,  nitric  oxide,  nitrogen 
dioxide,  carbon  monoxide,  and  total  oxidants. 
Definition  of  terms 

The  following  definitions  are  useful  in  the  proper  interpretation  of  the  data. 

1)  Hourly  average:  The  average  concentration  for  a  60-minute  period,  usually  the  clock 
hour.  This  is  an  integrated  value,  usually  read  from  the  strip  chart  by  eye.  Some  parts  of 
the  hour  have  lower  values  than  the  average,  some  higher.  There  are  24  such  values  for  each 
day  starting  with  hour  00  and  ending  with  hour  23. 

2)  Daily  maximum:  The  highest  concentration  reached  during  24  hours  of  sampling  and 
may  be  of  only  a  few  minutes  duration. 

3)  Daily  maximum  hourly  average:   This  is  the  highest  of  the  24  hourly  averages  re- 
ported during  the  day.   The  monthly  maximum  hourly  average  is  the  highest  such  value  reported 
during  the  month . 

4)  Daily  concentration:  The  recorded  value  is  the  reading  for  the  clock  hour  having 
the  highest  average  value.  The  clock  hour  is  any  one  of  24  60-minute  intervals  beginning 
and  ending  on  the  hour. 

Measurement  techniques 

1)  Total  oxidant  -  Refers  to  all  substances  which  may  be  present  in  the  atmosphere  that 
oxidize  potassium  iodide  (KI)  to  form  iodine  (I?),  such  as  ozone,  nitrogen  dioxide,  and 
organic  peroxides.   Specific  spectroscopic  measurements  indicate  ozone  is  the  principal  ox- 
idant in  California  air.   Concentrations  are  measured  with  continuous  recording  analyzers 
using  potassium  iodide  and  are  reported  as  parts  of  oxidant  per  million  parts  of  air  by 
colume.   In  Sacramento  the  instrument  is  colorimetric  and  is  standardized  against  ozone. 
Oxidant  measurements  for  Sacramento  have  been  printed  in  The  Clean  Air  Quarterly  for  each 
day  beginning  on  August  20,  1957. 

Note:   The  phenolphthaMn  method  for  measuring  oxidant  concentration  was  used  in  Sac- 
ramento until  August  24,  1960  at  which  time  the  potassium  iodide  method  came  into  use.   The 
phenolphthalein  method  was  as  follows:   Concentrations  were  determined  by  hand  analysis  of 
individual  samples  and  reported  as  parts  per  million  of  air  by  volume.   The  method  was 
standardized  against  hydrogen  peroxide.   Grab  samples  were  taken  once  a  day,  usually  between 
1000  and  1400  hours. 

2)  Nitrogen  oxides  (NO,  NO_,  NO  )  -  Continuous  measurement  of  NO~  is  performed  by  scrub- 
bing the  sample  air  with  an  absorbing  solution  (Saltzman  type  reagent)  in  an  automatic  re- 
cording analyzer.   Reaction  with  the  NO-  forms  a  red  AZO  dye  which  is  photometrically  measured 
and  recorded  on  the  strip  chart  as  ppm  N0_. 

NO  is  determined  by  first  passing  the  air  sample  through  an  acid  permanganate  bubble 

(oxidizer)  which  converts  NO  to  N0».   The  NO  plus  N0_  originally  present  is  recorded  as 

Oxides  of  Nitrogen  (NO  )  in  parts  per  million  of  air  by  volume.   The  NO  concentration  is 

taken  as  the  difference  between  the  NO  and  NO_  concentrations. 

x       2 

Daily  values  for  the  concentrations  of  NO  and  NO-  are  reported  in  the  Clean  Air  Quarterly 
beginning  on  October  4,  1960.   The  addition  of  Oxides  of  Nitrogen  to  the  published  data 

78 


tables  came  on  January  1,  1963.   Since  this  date  all  three  entries  have  been  made. 

3)  Carbon  monoxide  -  CO  is  measured  continuously  by  measuring  the  absorption  of  nondis- 
persed  infrared  radiation  (NDIR)  by  the  CO  nolecule.   Irradiation  from  a  hot  filament  is 
directed  to  a  detector  through  two  cells:  a  reference  cell  filled  with  dry  air  or  nitrogen, 
and  a  sample  cell  through  which  the  air  is  continuously  drawn.   The  difference  in  infrared 
absorption  between  the  reference  and  sample  cells  creates  a  pressure  difference  in  the  de- 
tector.  This  pressure  difference  is  transduced  to  CO  concentration  in  parts  per  million. 
The  principal  interferences  for  this  method  are  caused  by  water  vapor  and  carbon  dioxide. 
The  analyzer  is  of  the  recording  nondispersive  infrared  (Luft)  gype. 

4)  Hydrocarbons  -  Hydrocarbons  are  measured  with  continuous  flame  ionization  analyzers 
standardized  against  propane.   Results  are  reported  as  parts  of  methane  per  million  parts 
of  air. 

5)  Suspended  particulate  matter  -  Particulate  matter  samples  are  taken  by  automatic 
AISI  samplers  using  two-hour  sampling  cycles.   Values  are  calculated  by  comparing  light 
transmission  through  sampled  areas  and  blank  tape  and  are  expressed  as  COH  units*  per  1000 
linear  feet  of  air  passing  through  the  filter.   COH  units  less  than  1.0  indicate  relatively 
clean  air. 

Note:   Particulate  concentrations  for  dates  up  through  February  1967  are  tabulated  when 
readings  are  available  for  at  least  nine  (9)  2-hour  cycles.   Beginning  in  March  1967,  daily 
averages  are  tabulated  regardless  of  the  number  of  readings  reported. 
Dates  for  which  data  are  available 


1)  Oxidant 

a)  08/01/57  -  12/31/60 

b)  01/01/61  -  12/31/62 

c)  01/01/63  -  12/31/69 


Listed  as  daily  maximum  concentration 

Listed  as  daily  concentration 

Listed  as  maximum  hourly  average  concentration  by  day 

2)  Nitric  oxide  and  nitrogen  dioxide 

a)  10/01/60  -  12/31/60:   Listed  as  daily  maximum  concentration 

b)  01/01/61  -  12/31/62:   Listed  as  daily  concentration 

c)  01/01/63  -  12/31/69:   Listed  as  maximum  hourly  average  concentration  by  day 

3)  Oxides  of  nitrogen 

01/01/63  -  12/31/69:   Listed  as  maximum  hourly  average  concentration  by  day 

4)  Hydrocarbons 

01/01/65  -  12/31/69:   Listed  as  maximum  hourly  average  concentration  by  day 

5)  Carbon  monoxide 

Listed  as  daily  maximum  concentration 
No  data  reported 

Listed  as  daily  concentration 

Listed  as  maximum  hourly  average  concentration  by  day 


a)  01/01/60  -  06/30/60 

b)  07/01/60  -  12/31/61 

c)  01/01/62  -  12/31/62 

d)  01/01/63  -  12/31/69 


6)  Particulate  matter 

04/01/58  -  12/31/69:   Listed  as  daily  mean  COH  values  by  the  AISI  method 
Note:   Particulate  concentrations  for  dates  04/01/58  -  02/28/67  are  tabulated 


*  Note:   One  "COH  Unit"  is  the  quantity  of  light  scattering  solids  which  is  capable  of 
producing  an  optical  density  of  0.01  for  light  at  0.4y. 
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only  when  readings  for  nine  (9)  2-hour  cycles  are  available.   Beginning  in  March 
1967,  daily  averages  are  tabulated  regardless  of  the  number  of  readings  reported. 
Location  of  sampling  stations 

Most  of  the  pollutant  sampling  which  is  reflected  in  our  data  has  been  measured  at  the 
downtown  station  (1230  J  Street)  .   Oxidants  have  been  continuously  measured  there  since 
September  16,  1961;  particulate  matter  since  June  8,  1962;  NO  and  N0_  since  September  16, 
1961;  and  carbon  monoxide  since  October  5,  1961.   In  1965  hydrocarbons  were  added.   The  height 
of  the  sampling  point  is  about  25  feet  above  the  ground. 

The  data  which  are  on  punched  cards  begin  with  January  1,  1961.   Thus,  any  pollution 
data  for  the  period  January  1,  1961  -  September  15,  1961  was  measured  at  2221  Stockton  Blvd. 
where  oxidants  and  particulates  were  sampled  during  the  early  years  of  air  contaminate  mon- 
itoring.  The  sampling  point  at  the  latter  location  was  approximately  3  feet  above  ground 
level. 
Critique  of  instrumentation  used  to  obtain  pollution  data* 

(1)  Oxidant.   The  continuously  recording  coiorimetric  analyzers  used  to  measure  total 
oxidant  have  been  calibrated  semi-annually  with  ozone-air  mixtures  so  that  the  oxidant  con- 
centrations in  the  air  are  registered  as  apparent  0_  concentrations.   Over  a  4-year  period, 
67%  of  the  instruments  in  use  in  California  have  been  found  with  -10%  and  83%  within  -20% 
of  the  previous  calibration. 

The  principal  oxidant  species  besides  0  which  oxidize  potassium  iodide  to  I_  are  NO-, 
PAN,  organiz  peroxides,  Cl_,  H_02,  etc.   The  sensitivity  of  the  KI  method  to  interference  by 
N0„  increases  with  the  concentration  of  the  KI  solution,  and  with  a  20%KI  solution  (the 
maximum  concentration  used)  30  to  40%  of  the  N0~  present  is  measured  as  ozone.   PAN  concen- 
trations in  the  atmosphere  can  be  as  large  as  .04  ppm  and  50%  of  the  PANs  are  recorded  as 
apparent  0,  concentrations.   Hence  these  species  cause  interference  as  well  as  doubts  about 
the  species  content  in  the  group  of  oxidants.   S0?  is  found  to  interfere  negatively,  result- 
ing in  a  decrease  in  the  oxidant  reading  by  an  amount  equal  to  the  concentration  of  S0? 
present  in  the  sample.   In  general,  the  performance  of  these  and  other  instruments  based  on 
wet  chemical  analysis  appears  to  depend  greatly  on  the  diligence  and  skill  of  the  person 
tending  them.   Reagent  quality  and  concentrations,  air  flow,  and  air-solution  contact  are 
regulated  by  non-reliable  automated  wet-channel  devices  which  require  close  attention.   The 
instrument  in  use  in  Sacramento  is  made  by  Beckman  instruments. 

(2)  Oxides  of  nitrogen  are  now  measured  by  a  continuous  automatic  coiorimetric  analyzer 
(Beckman).   The  Beckman  instruments  used  in  California  have  been  calibrated  semi-annually 
for  N0».   In  stability  tests  67%  of  the  analyzers  were  found  to  be  within  -10%  of  the  pre- 
vious calibration  and  85%  with  -20%.   The  lower  limit  for  meaningful  concentration  measure- 
ments was  0.01  ppm.   Calibration  for  NO  has  been  incomplete  because  of  the  lock  of  a  re- 
liable NO  source.   Therefore,  all  measurement  data  actually  represent  N0?  equivalent  con- 
centrations. 

In  the  parallel  mode  (used  in  Sacramento)  NO-  is  lost  in  the  oxidizer  and  only  about  80% 
of  the  NO  present  is  converted  to  N0».   Also  50%  of  the  concentration  of  PANs  is  measured 


kThis  section  is  based  on  Project  Clean  Air,  University  of  California,  Task  Force  Assess- 
ments, Volume  4. 
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as  N0_.   Strict  adherence  to  the  Saltzmann  procedure  is  required  to  avoid  important  inter- 
ferences by  SO^  and  NO  with  the  N0?  determination  in  air. 

Furthermore,  no  information  has  appeared  about  the  reliability  of  the  NO  ■*  NO     conver- 
sion in  KM  ) ,  acid  solution,  nor  of  the  interferences  with  this  process, 
n  4  r 

(3)  Carbon  monoxide  is  presently  measured  by  the  Beckman  Non-dispersive  Infrared  Anal- 
yzer (NDIR) .   This  device  is  believed  to  be  relatively  accurate  and  reliable.   However,  if 
proper  filtering  techniques  are  not  used  in  the  analyzers,  H_0  and  C0?  cause  the  principal 
interferences  and  to  a  lesser  degree  CH  and  C?H  .   The  major  drawback  with  NDIR  analysis  is 
the  rather  ambiguous  air  sampling  technique.   Air  samples  are  collected  in  aluminized  bags 
and  then  analyzed.   A  more  direct  means  of  air  intake  to  the  analyzing  instrument  at  a  con- 
stant regulated  flow  rate  would  be  preferable.   The  lower  limit  for  meaningful  concentration 
measurements  is  about  1  ppm. 

(4)  Total  hydrocarbon  is  now  measured  by  Hydrogen  flame  ionization  analysis.   No  per- 
formance data  seem  to  be  available  for  appraisal  of  the  value  of  the  instruments  presently 

in  use.   However,  once  proper  flame  conditions  and  precise  maintenance  of  flow  rate  have  been 
established,  the  breakdown  problems  seem  minimal.   The  lower  limit  of  meaningful  concentration 
measurement  seems  to  be  about  1  ppm  of  methane  equivalent. 

(5)  Particulates  -  All  the  particulate  concentrations  reflected  in  the  Sacramento  data 

are  recorded  in  COH  units  per  1000  linear  feet  of  air  sampled.   While  these  measurements  can 

3 
be  related  approximately  to  the  concentrations  in  ug/m  of  particulates  in  the  air  sample, 

they  depend  strongly  on  the  number,  size,  color,  and,  in  particular,  on  the  phase  (liquid  or 
solid)  of  the  suspended  matter  as  well  as  on  the  grade  of  filter  material  used. 
Meteorological  data 

Meteorological  data  were  taken  from  the  monthly  issues  of  the  U.  S.  Department  of  Com- 
merce, Weather  Bureau,  in  the  form  of  the  Local  Climatological  Data  (Supplement)  January 
1961  -  December  1969.   Following  is  a  summary  of  the  dates  for  which  the  various  data  are 
available  and  an  explanation  of  reporting  units. 

1)  Average  dry-bulb  temperature  (°F) 

January  1,  1961  -  December  31,  1969:  data  available 

T    +  T 
Average  appears  to  be  max    min  for  each  day. 

2 

2)  Precipitation  (water  equivalent)  in  inches 

January  1,  1961  -  March  31,  1962:  no  data 

April  1,  1962  -  April  30,  1962:  data  available 

May  1,  1962  -  August  31,  1962:  no  data 

September  1,  1962  -  September  30,  1962:  data  available 

October  1,  1962  -  September  30,  1962:  no  data 

October  1,  1963  -  November  30,  1963:  data  available 

December  1,  1963  -  December  31,  1963:  no  data 

January  1,  1964  -  May  31,  1964:  data  available 

June  1,  1964  -  August  31,  1964:  no  data 

September  1,  1964  -  December  31,  1969:  data  available 

3)  Resultant  wind  direction  measured  from  north  taken  clockwise.   (Northerly  =  360, 

Easterly  =  180,  Westerly  270,  0  =  calm). 

January  1,  1961  -  December  31,  1969:  data  available 
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4)  Resultant  wind  speed  (mph) .   Resultant  wind  movement  is  the  vectorial  sum  of  all 

surface  wind  directions  and  speeds  for  each  day.   It  is  obtained  by  resolving  each 
hourly  wind  observation  into  components  from  north  and  east,  summing  over  the  day, 
and  reconverting  the  summed  components  into  a  single  vector.   The  daily  resultant 
wind  speed  is  obtained  by  dividing  the  total  movement  by  24. 

All  meteorological  data  reported  were  measured  at  the  Sacramento  Municipal  Airport.  (Name 

changed  October  1,  1969  to  Executive  Airport). 

Statistical  analysis  of  the  data 

After  the  foregoing  data  were  transferred  to  punched  cards,  a  careful  error  check  was 

made,  followed  by  an  extensive  study  of  the  elementary  statistics  which  could  be  formed  from 

these  data  including  the  arithmetic  mean,  the  standard  deviation,  maximum  and  minimum  values, 

and  correlation  coefficients  between  all  input  parameters.   These  calculations  were  made  for 

the  following  groupings  of  data: 

(1)  Each  individual  month  (108  months  -  no  correlations) 

(2)  Each  individual  year  (1961  to  1969) 

(3)  All  nine  years  together 

(4)  All  nine  years  together,  Sunday  only 

(5)  January  -  March 

(6)  January  -  March,  with  wind  restrictions* (any  day  with  wind  directions  between 

east  and  north  were  excluded) 

(7)  April  -  June 

(8)  April  -  June,  with  wind  restrictions* 

(9)  July  -  September 

(10)  July  -  September,  with  wind  restrictions* 

(11)  October  -  December 

(12)  October  -  December,  with  wind  restrictions* 

(13)  December  -  February 

(14)  March  -  May 

(15)  June  -  August 

(16)  September  -  November 

A  few  words  regarding  the  limitations  of  these  data  are  in  order.   The  pollution  data 
were  available  in  the  form  of  daily  values  of  maximum  hourly  averages.   For  some  applications 
relating  to  public  health  and  human  comfort,  such  information  is  quite  useful.   On  the  other 
hand,  these  data  severely  limits  our  ability  to  make  detailed  verification  tests  for  model 
development.   For  this  application,  average  daily  values  would  be  preferable  and  direct  access 
to  hourly  averages  would  be  ideal.   However,  these  data  are  still  useful  as  numerical  base- 
line information  in  spite  of  the  limitations.   In  this  section,  we  shall  give  the  principal 
results  of  our  analysis.   We  have  not  yet  subjected  these  data  to  more  advanced  statistical 
analyses,  such  as  the  higher  moments,  probability  and  joint  probability  density  distribution, 
or  spectral  and  cross-spectral  analyses.   In  view  of  the  fact  that,  in  general,  a  daily  hourly 
maximum  for  any  one  pollutant  will  not  occur  at  the  same  hour  as  another,  any  correlation- 
type  analysis  must  be  viewed  with  caution.   It  should  also  be  remembered  that  the  particulate 
data,  unlike  the  rest,  exist  in  the  form  of  daily  mean  hourly  averages. 
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TABLE  2-1. 

MEAN  VALUES  1961 

-  1969 

Quantity 

Unit 

Average 

Maximum 

Standard 

Value 

Value 

Deviation 

Oxidant 

ppm 

0.055 

0.290 

0.0645 

NO 

ppm 

0.113 

1.50 

0.145 

N02 

ppm 

0.061 

0.50 

0.042 

NO 

X 

ppm 

0.158 

1.19 

0.159 

Hydrocarbon 

ppm 

6.460 

33.0 

4.326 

Particulates 

COH 

0.506 

2.80 

0.378 

CO 

ppm 

6.700 

55.0 

6.472 

Temperature 

°F 

60.2 

91.0 

12.2 

Wind  speed 

mph 

5.90 

22.0 

3.96 

Wind  direction 

deg. 

222 

--- 

75.0 

Note  that  temperature  and  wind  velocity  are  also  averaged  over  the  entire  day. 

Inspection  of  Table  2-1  shows  that  the  maximum  hourly  pollutant  concentration  is  a 
highly  variable  statistic  since  most  of  the  standard  deviations  are  of  the  same  order  of 
magnitude  as  the  mean  value. 

Sunday-only  data,  1961  -  1969  -  In  order  to  begin  to  investigate  the  relationship  be- 
tween human  activities  and  pollution  concentration,  we  computed  the  mean  values  for  all 
Sunday  data  over  the  same  period.   The  results  of  this  analysis  are  given  below. 


TABLE  2-2.   POLLUTANT  CONCENTRATIONS  FOR  SUNDAYS  ONLY,  1961  -  1969 


Pollutant 

Uni 

Oxidant 

ppm 

NO 

ppm 

N02 

ppm 

NO 

X 

ppm 

Hydrocarbons 

ppm 

Particulates 

COH 

CO 

ppm 

Mean 

Standard 

Percent 

Value 

Deviation 

Reduction 

0.0497 

0.0282 

9.1 

0.0845 

0.123 

25.7 

0.0429 

0.0317 

29.5 

0.125 

0.150 

20.9 

6.040 

4.042 

6.5 

0.373 

0.321 

26.3 

5.624 

5.979 

16.1 

All  the  primary  pollutants  show  reductions  compared  to  the  mean  values  computed  for  all 
days.   The  largest  reductions  appear  for  N0?,  particulates,  NO,  and  NO  with  lesser  reduc- 
tions for  the  remainder.   This  result  implies  that  those  few  pollutants  are  those  which  are 
most  directly  related  to  human  activities  on  a  24-hour  time  scale. 

Monthly  variation  -  A  tabulation  of  the  mean  monthly  values  of  the  various  parameters 
is  given  below. 
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TABLE  2-3 

MONTHLY  AVERAGES 

1961 

-  1969 

Month 

Oxidant 

NO 

N02 

NO 

X 

HC 

Part 

CO 

Temp. 

Wind 
Speed 

Wind 
Direction 

JAN. 

0.0358 

0.163 

0.0645 

0.203 

7.306 

0 

.762 

7.983 

43.65 

4.90 

206.6 

FEB. 

.0364 

.1304 

.0586 

.1813 

7.435 

.5853 

7.264 

49.70 

5.30 

218.7 

MAR. 

.0425 

.0740 

.0526 

.1211 

6.404 

.4308 

5.720 

52.63 

6.14 

215.8 

APRIL 

.0471 

.0421 

.0438 

.0830 

5.104 

.3002 

4.556 

57.84 

6.91 

230.0 

MAY 

.0554 

.0337 

.0448 

.0737 

4.435 

.2947 

3.451 

63.92 

7.21 

228.9 

JUNE 

.0563 

.0336 

.0488 

.0718 

3.971 

.2740 

3.449 

70.49 

7.45 

227.4 

JULY 

.0728 

.0299 

.0572 

.0792 

4.000 

.2984 

3.711 

74.87 

6.93 

213.3 

AUG. 

.0687 

.0441 

.0542 

.0902 

5.441 

.3557 

4.472 

75.14 

6.11 

218.4 

SEPT. 

.0711 

.0643 

.0647 

.1181 

5.440 

.5020 

5.449 

71.45 

6.09 

235.2 

OCT. 

.0605 

.2202 

.1000 

.3240 

10.72 

.7902 

13.81 

63.86 

4.85 

239.2 

NOV. 

.0350 

.2357 

.0768 

.2881 

9.236 

.7743 

11.64 

53.31 

4.55 

218.4 

DEC. 

.0240 

.2075 

.0637 

.2363 

7.094 

.6991 

8.307 

44.84 

4.61 

210.4 

These  numbers  are  plotted  in  Figures  2-1  and  2-2.   As  can  be  seen,  the  average  monthly 
pollution  concentrations  fall  into  two  patterns.   The  majority  of  these  pollutants  show  a 
maximum  in  October  (November  for  NO)  and  a  minimum  in  June  (July  for  NO,  April  for  NO~)  •   This 
pattern  closely  resembles  the  reciprocal  of  wind  speed  which  shows  a  maximum  in  November 
and  minimum  in  June  (Figure  2-4)  so  that  one  would  conclude  that  dilution  due  to  the  action 
of  the  wind  is  a  dominant  process  for  these  pollutants.  Oxidant  displays  an  entirely  dif- 
ferent pattern.   The  annual  maximum  is  reached  in  July,  followed  by  decreasing  values  during 
the  period  when  the  other  pollutants  are  increasing.   The  minimum  for  oxidant  is  reached  in 
December  close  to  the  time  of  the  maximum  for  the  other  pollutants.   This  different  be- 
havior on  the. part  of  oxidant  most  certainly  is  a  reflection  of  the  importance  of  its  photo- 
chemical origin.   The  generation  of  oxidant  is  a  function  of  the  intensity  of  solar  radia- 
tion and  the  concentration  of  the  reagents  which  participate  in  this  complex  reaction  (NO, 
N0„,  hydrocarbons).   Evidently,  conditions  are  optimum  in  July  with  decreasing  generation  of 
oxidant  thereafter  in  spite  of  increasing  concentrations  of  the  input  reagents  due  to  de- 
creasing radiation.   The  remaining  meteorological  parameters  are  plotted  in  Figure  2-3. 

Seasonal  statistics  -  Since  air  pollution  is  more  of  a  seasonal  phenomena  than  a 
monthly  one,  we  computed  the  same  statistics  as  before  on  a  seasonal  basis.   We  chose  the 
periods  January  -  March,  April  -  June,  July  -  September,  and  October  -  December  for  this 
analysis.   Below  are  the  results.   These  values  show  essentially  the  same  pattern  as  our 
previous  analysis  and  are  not  plotted.  (See  Table  2-4) 

Correlation  analysis  -  We  have  also  subjected  these  data  to  a  correlation  analysis.   In 
a  sense,  the  limitations  of  daily-maximum  hour-concentrations  are  an  advantage  here.   The 
rate  constants  in  the  complex  photochemical  process  are  such  that  peak  values  in  the  various 
participants  occur  at  different  times.   However,  a  bad  day  for  one  pollutant  would  be  ex- 
pected to  be  bad  for  one  of  the  other  chemical  actors  but  with  the  maximum  concentration 
occurring  at  a  different  time.   Hence,  correlations  between  the  various  pollutants  without 
regard  for  the  time  of  their  maximum  values  may  be  more  meaningful  than  if  the  analysis  were 
done  hour-by-hour  without  consideration  of  time-lags  among  the  pollutants.   Table  2-5  shows 
the  correlation  matrix  for  all  data  available  for  1961  -  1969. 
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Figure  2-1.   Monthly  averages  of  the  daily  maximum 
concentration  of  Total  Oxidant,  Nitric 
Oxide,  Nitrogen  Dioxide  and  Oxides  of 
Nitrogen.   1961-1969. 
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Figure  2-2.   Monthly  averages  of  the  daily  maximum 

concentration  of  Carbon  Monoxide,  Hydro- 
carbons and  Particulates.   1961-1969. 
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Figure  2-4.   Monthly  average  of  the  reciprocal  of  wind 
speed.   1961-1969. 
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TABLE  2-4. 

SEASONAL 

STATISTICS  1961-1969 

Oxidant 

Jan. 

Mean 

-Mar. 

0 

Apri 
Mean 

1-June 

0 

Jul 

Mear 

y-Sept. 

L         0 

Oct.- 

Mean 

Dec. 

0 

Oxidant 

.0401 

.0241 

.0536 

.0264 

.072* 

1     .0360 

.0399 

.0319 

NO 

.1256 

.1173 

.0392 

.0359 

.0485 

1     .0558 

.2220 

.2023 

N02 

.0565 

.0346 

.0471 

.0226 

.059] 

.0396 

.0805 

.0540 

NO 

X 

.1704 

.1292 

.0766 

.0424 

.0978 

!     .0737 

.2856 

.2234 

HC 

7 

.082 

3.739 

4.538 

2.302 

5 

.011 

2.559    8 

.058 

5 

.821 

Particulate 

.5885 

.3832 

.2879 

.1506 

.386J 

1     .2379 

.7518 

.4637 

CO 

7 

.035 

4.916 

3.799 

2.225 

4 

.586 

3.554    11 

.290 

9 

.527 

Temp. 

48 

.70 

6.26 

64.01 

8.25 

73 

.85 

5.39    54 

.01 

9 

.59 

Speed 

5 

.40 

4.18 

7.16 

3.71 

6, 

.38 

3.33     4 

.67 

4 

.15 

Direction 

215 

.3 

88.9 

229.3 

51.37 

222. 

.3 

43.5    222 

.9 

99 

.3 

TABLE  2-5.   CORRELATION  COEFFICIENTS  1961-1969 

Ox 

idant 

NO 

N02 

NO 

X 

HC     1 

'art, 

CO     Temp . 

Speed 

Direction 

Oxidant 

1 

-.035 

.286 

.036 

.104 

.117 

108     .534 

-.126 

-.110 

NO 

.035 

1 

.514 

.981 

.806 

.690 

812    -.304 

-.409 

.070 

N02 

.286 

.514 

1 

.604 

.836 

.571 

847     .018 

-.375 

.002 

NO 

X 

.036 

.981 

.604 

1 

.836 

.741 

847    -.246 

-.442 

.081 

HC 

.104 

.806 

.836 

.836 

1 

.781 

866   -.144 

-.480 

.086 

Part. 

.117 

.690 

.571 

.741 

.781 

1 

754    -.257 

-.551 

.051 

CO 

.108 

.812 

.847 

.847 

.866 

.754 

1     -.157 

-.449 

.067 

Temp. 

.534 

-.304 

.018 

-.246 

.144 

.257 

-. 

157      1 

.120 

.076 

Speed 

.126 

-.409 

-.375 

-.442 

.480 

.551 

-. 

449     .120 

1 

.067 

Direction- 

.110 

.070 

.002 

.081 

.086 

.051 

067     .076 

.067 

1 

Due  to  the  large  number  of  degrees  of  freedom,  the  great  majority  of  these  correlations  are 
significant  at  the  1%  level  or  better.   The  following  features  are  particularly  noteworthy. 
For  oxidant,  the  correlation  coefficients  are  unifornly  low  except  for  temperature.   Among 
the  remaining  pollutants,  the  correlations  are,  without  exception,  high.   The  highest  cor- 
relations are  between  NO  and  NO  (0.981),  CO  and  hydrocarbon  (0.866),  CO  and  NO.,  (0.847),  and 

X  z 

CO  and  NO   (0.847).   Between  pollutants  and  meteorological  perameters,  only  wind  speed  shows 
a  consistent  pattern.   With  the  exception  of  oxidant,  the  correlation  between  wind  speed  and 
pollution  concentration  tends  to  be  about  -0.5,  with  particulates  yielding  the  largest 
negative  correlation  (-0.551).   Oxidant  is  the  only  pollutant  that  shows  a  significant  cor- 
relation with  temperature  and  no  significant  correlations  are  found  with  wind  direction. 

Most  of  these  results  are  unsurprising.   The  primary  source  for  all  these  pollutants  is 
the  automobile,  so  they  should  be  correlated.   Similarly,  the  diluting  effect  of  the  wind 
would  be  expected  to  be  reflected  as  a  negative  correlation.   The  low  oxidant  correlations 
are  surprising  at  first  glance  because  it  is  well  known  from  the  Los  Angeles  experience  that 
high  oxidant  hays  tend  to  be  high  in  the  other  major  pollutants.   One  hypothesis  to  explain 
the  Sacramento  data  is  that  during  the  late  fall  months  when  most  of  the  pollutants  are 
building  up  due  to  lessening  wind  speed,  oxidant  decreases  due  to  decreasing  radiation.   This 
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effect  is  especially  pronounced  in  Sacramento  because  of  a  maximum  in  fog  in  the  late  fall 
and  early  winter.   However,  sunny  days  do  occur,  even  during  the  peak  of  the  fog  season.   On 
these  days  oxidant  could  be  high,  along  with  temperature.   Hence,  the  correlation  between 
temperature  and  oxidant  should  persist  throughout  the  year. 

To  test  this  hypothesis,  we  performed  the  correlation  analysis  for  the  various  seasons, 
as  defined  above.   For  oxidant,  we  obtained  the  following  results. 

TABLE  2-6.   SEASONAL  CORRELATION  COEFFICIENTS  WITH  OXIDANT  1961-1969 
Quantity Jan.  -Mar. April -June July-Sept .       Oct. -Dec. 


NO 

0.1399 

-0.0301 

0.1515 

0.2516 

N02 

0.3438 

0.2220 

0.3422 

0.5336 

NO 
x 

0.1760 

0.0557 

0.2331 

0.3631 

HC 

0.2284 

0.1292 

0.2670 

0.5020 

Part. 

0.3229 

0.2889 

0.3115 

0.4160 

CO 

0.3464 

0.1668 

0.2197 

0.4121 

Temp. 

0.2514 

0.3833 

0.2970 

0.5525 

Speed 

-0.1509 

-0.1997 

-0.3384 

-0.1435 

Direction 

0.0233 

-0.0725 

-0.0632 

0.0219 

As  can  be  seen  by  comparing  these  values  with  those  in  the  previous  table,  the  correlation  be- 
tween oxidant  and  the  other  major  pollutants  is  higher  in  almost  every  case  than  for  the  pre- 
vious annual  analysis.   The  correlations  during  the  late  fall  are  particularly  high.   For  the 
other  pollutants  and  meteorological  parameters,  the  effect  of  seasonal  grouping  is  not  as 
pronounced.   For  instance,  the  following  results  were  obtained  for  carbon  monoxide. 


TABLE  2-7.   SEASONAL  CORRELATION  COEFFICIENTS  WITH  CARBON  MONOXIDE  1961-1969 
Quantity Jan.  -Mar. April -June July-Sept . Oct.  -Dec  . 


Oxidant 

.3464 

NO 

.7235 

N02 

.4641 

NO 
x 

.7602 

HC 

.8205 

Part. 

.6947 

Temp. 

-.0938 

Speed 

-.4743 

Direction 

.0114 

.1668 
.5695 
.3702 
.6558 
.6012 
.5126 
,0725 
.3796 
.1687 


.2197 
.6823 
.4191 
.7349 
.8260 
.6799 
.1109 
.4394 
.3304 


.4121 
.7933 
.5859 
.8369 
.8888 
.7477 
.2352 
.4321 
.0844 
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For  particulate,  the  seasonal  analysis  yielded  the  following  values. 

TABLE  2-8.   SEASONAL  CORRELATION  COEFFICIENTS  WITH  PARTICULATES  1961-1969 
Quantity Jan.  -Mar. April -June July-Sept. Oct .  -Dec  , 


Oxidant 

.3229 

NO 

.5414 

N02 

.4909 

NO 

X 

.5707 

HC 

.6462 

CO 

.6947 

Temp. 

-.3010 

Speed 

-.5567 

Direction 

-.0114 

.2889 
.3361 
.3953 
.4742 
.4024 
.5126 
.2351 
.4494 
.1524 


.3115 
.5658 
.4927 
.6810 
.7221 
.6799 
.1708 
.5078 
.3348 


.4160 
.6703 
.6230 
.7604 
.8102 
.7477 
.0758 
.5639 
.0779 


Finally,  the  seasonal  correlation  coefficients  for  wind  speed  are  the  following: 

TABLE  2-9.   SEASONAL  CORRELATION  COEFFICIENTS  WITH  WIND  SPEED  1961-1969 
Quantity Jan.  -Mar. April -June July-Sept. Oct .  -Dec  , 


Oxidant 

-0.1509 

NO 

-0.3660 

N02 

-0.3017 

NO 

X 

-0.3802 

HC 

-0.4683 

Part. 

-0.5567 

CO 

-0.4743 

Temp. 

0.1707 

Direction 

0.0634 

Temporal  variation 

-  As  we  1 

-0.1997 
-0.3730 
-0.4081 
-0.4729 
-0.3304 
-0.4494 
-0.3796 
-0.1337 
-0.0619 


-0.3384 
-0.3817 
-0.3891 
-0.4745 
-0.5209 
-0.5078 
-0.4394 
-0.3330 
-0.2113 


-0.1435 
-0.4141 
-0.3378 
-0.4735 
-0.4504 
-0.5639 
-0.4321 
0.1047 
0.0907 


examining  these  statistics  is  to  study  whatever  temporal  trends  may  exist  in  air  quality. 
We  have  found  it  useful  in  attempting  to  assess  the  existence  of  trends  to  keep  in  mind  the 
major  processes  which  effect  air  quality.   During  the  decade  1960-70,  in  which  population 
increased  by  22%  and  energy  consumption  by  50%,  the  number  of  motor  vehicals  increased  by 
approximately  30%.   In  the  absence  of  emission  control  efforts,  and  all  other  factors  being 
the  same,  one  would  expect  air  pollution  concentrations  to  increase  by  about  3  or  4%  per 
year  during  the  period  in  question.   We  shall  presently  examine  the  variability  of  meteoro- 
logical parameters.   The  control  effort  for  automobile  emission  began  about  1962  when  crank- 
case  systems  began  to  be  installed.   The  mechanism  is  thought  to  have  become  a  significant 
factor  in  air  quality  by  about  1965.   This  control  approach  is  effective  primarily  for 
hydrocarbons  so  that  impact  would  be  expected  only  for  the  photochemical  pollutants. 

In  the  following  sections  we  present  the  elementary  temporal  statistics.   We  have  not 
subjected  the  data  to  an  elaborate  regression  analysis.   As  will  become  apparent  in  the  dis- 
cussion, we  do  not  feel  that  the  evidence  warrants  such  a  treatment.   Table  2-10  shows  an- 
nual averages  for  the  pollution  and  meteorological  statistics  and  their  standard  deviations 
for  each  year.   The  standard  deviations,  it  will  be  noted,  are  of  similar  magnitude  to  the 


91 


averages  in  all  cases  for  the  pollutants.  In  the  series  of  figures  which  follow,  the  annual 
average  pollutant  concentrations  are  shown  plotted  in  groups  with  parameters  which  have  pre- 
viously been  shown  to  have  significant  correlations. 

TABLE  2-10.   ANNUAL  AVERAGES  OF  POLLUTION  CONCENTRATIONS  AND  METEOROLOGICAL 
PARAMETERS  1961-1969 
1961     1962     1963     1964     1965     1966     1967     1968     1969 


Oxidant 

.0371 

.0553 

.0457 

.0597 

.0620 

.0494 

.0528 

.0423 

.0561 

0 

.0310 

.0332 

.0284 

.0364 

.0341 

.0310 

.0344 

.0271 

.0332 

NO 

.1543 

.0866 

.1163 

.1073 

.1135 

.1157 

.1250 

.0960 

.1019 

a 

.2153 

.1045 

.1423 

.1330 

.1514 

.1317 

.1496 

.1267 

.1281 

N02 

.0697 

.0655 

.0724 

.0625 

.0672 

.0504 

.0613 

.0464 

.0570 

0 

.0559 

.0446 

.0444 

.0431 

.0402 

.0391 

.0400 

.0261 

.0306 

NO 
x 

.1699 

.1540 

.1681 

.1587 

.1731 

.1356 

.1454 

0 

.1573 

.1568 

.1731 

.1596 

.1748 

.1401 

.1444 

HC 

7.396 

6.137 

7.045 

6.129 

5.853 

o 

4.734 

3.927 

5.156 

3.782 

3.939 

Part. 

.6063 

.5268 

.5008 

.4874 

.5183 

.5112 

.5970 

.4635 

.3587 

a 

.3871 

.4191 

.4182 

.3670 

.3678 

.3324 

.432 

.3346 

.2681 

CO 

6.597 

6 

.177 

6 

.740 

7.311 

6.525 

7.268 

6.297 

6.662 

0 

6.349 

4 

.841 

6 

.587 

7.264 

6.139 

7.159 

5.974 

7.048 

Temp. 

59 

.69 

58.77 

58 

.40 

60 

.22 

59.92 

61.61 

60.95 

61.13 

61.01 

a 

12 

.88 

11.95 

12 

.09 

11 

.28 

11.72 

12.00 

13.39 

12.12 

12.41 

Speed 

5 

.76 

5.42 

4 

.60 

7 

.97 

6.14 

6.00 

5.26 

5.94 

6.05 

0 

4 

.04 

3.93 

3 

.58 

3 

.71 

4.12 

3.88 

3.71 

3.90 

3.98 

(Speed) 

1 

.1736 

.1845 

.2174 

.1255 

.1629 

.1667 

.1901 

.1683 

.1653 

Direct. 

221 

.7 

224.8 

224 

.4 

224 

.2 

218.8 

225.5 

2.6.8 

227.4 

218.4 

a 

76 

.2 

70.3 

70 

.0 

72 

.0 

81.4 

75.7 

76.5 

79.0 

75.03 

In  Figure  2-5  the  annual  data  are  shown  for  hydrocarbons,  N0_  and  inverse  wind  speed. 
The  data  period  for  hydrocarbon  is  too  short  to  permit  conclusions,  but  a  downward  trend  is 
suggested.   The  NO,,  data  likewise  show  an  apparent  downward  trend  throughout  the  entire  nine- 
year  period.   The  small  positive  correlation  with  inverse  wind  speed  can  easily  be  seen. 
Figure  2-6  shows  the  data  for  oxidant,  temperature  and  inverse  wind  speed.   Disappointingly, 
no  trend  can  be  discerned  for  oxidant,  and  negligible  interrelationships  can  be  detected  be- 
tween these  variables.   Figure  2-7  shows  the  same  analysis  for  carbon  monoxide,  particulate 
matter  and  inverse  wind  speed.   No  trend  is  evident  for  carbon  monoxide.   Particulates  do 
indicate  a  decrease  in  the  last  three  years  of  the  data  period,  but  this  can  be  seen  to  be 
highly  correlated  with  a  decrease  in  the  inverse  wind  speed  at  the  same  time. 

Our  previous  analysis  indicated  that  correlations  among  the  pollution  and  meteorological 
variables  were  more  apparent  if  the  analysis  was  done  on  a  seasonal  basis.   Accordingly,  we 
have  examined  temporal  variation  on  a  mean  monthly  basis.   Tables  2-11  and  2-12  show  the 
monthly  averages  for  June  and  October  for  each  year  in  the  study  period.   Figure  2-8  shows  a 
plot  of  the  June  data  for  oxidant,  temperature  and  inverse  wind  speed.   No  trend  is  evident 
for  oxidant  and  interrelationships  seem  nil.   The  previous  analysis  had  indicated  weak 
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positive  correlations  to  exist  among  these  variables  but  that  is  not  apparent  for  this  month. 
The  year  1965  is  particularly  puzzling  in  that  the  highest  June  oxidant  concentration  is 
reached  in  this  year  despite  the  lowest  values  of  temperature  and  inverse  wind  speed.   It 
could  be  argued  that  this  maximum  is  the  result  of  increasing  numbers  of  motor  vehicles  with, 
as  yet,  ineffective  emission  control.   The  decrease  after  1965  would  then  reflect  the  in- 
creased effectiveness  of  the  control  technology.   Figure  2-9  shows  the  June  data  for  carbon 
monoxide,  particulates  and  inverse  wind  speed.   No  trend  is  apparent  and  obvious  correlations 
can  be  seen  between  these  parameters. 

TABLE  2-11.   JUNE  MONTHLY  AVERAGES  1961-1969 

1961    1962    1963    1964    1965    1966    1967  1968  1969 

Oxidant    0.0460    .0625    .0282    .0589    .0934    .0604    .0662  .0460  .0487 

CO           --    3.759   2.700   4.103   3.233   4.133   3.200  3.172  3.310 

Part.       .4286    .2567    .2033    .2464    .2300    .3321    .2960  .2714  .2033 

Temp.     74.3    70.03   68.70   68.27   67.30   72.80   69.60  74.17  69.33 

Speed      7.06    7.24    6.78    7.07    9.61    6.61    6.92  7.04  8.79 

(Speed)"1    .1416    .1381    .1474    .1414    .1040    .1512    .1445  .1420  .1137 

TABLE  2-12.   OCTOBER  MONTHLY  AVERAGES  1961-1969 
1961    1962    1963    1964    1965    1966    1967    1968    1969 
Oxidant     .0321    .0493    .0568    .0865    .0800    .0694    .0766    .0559    .0397 

CO  --    12.82    9.129   15.25    19.06   13.21    17.43  13.40  10.16 

Part.       .9192    .5968    .6821    .9862   1.026    .7903    .9964   .7214    .3948 

Temp.     62.58   61.00   64.48   66.90   66.22   64.26   65.74  62.35  61.26 

Speed      6.26    4.87    4.12    5.70    3.75    5.21    4.20  4.39    5.21 
(Speed)"1    .1597    .2053    .2427    .1754    .2667    .1919    .2380    .2278    .1919 

The  next  series  of  figures  show  the  same  analysis  for  the  month  of  October.   Data  for 
oxidant,  temperature  and  inverse  wind  speed  are  given  in  Figure  2-10.   The  oxidant  ^profile 
for  Ocrober  can  be  seen  to  follow  a  pattern  very  similar  to  that  for  June  except  that  the 
maximum  is  reached  one  year  sooner.   On  the  face  of  it,  this  could  interpreted  as  an  en- 
couraging indication  that  the  emission  control  program  is  having  a  significant  effect  since 
there  is  little  meteorological  reason  for  events  in  October  to  be  correlated  with  those  in 
June,  on  the  average.   However,  further  examination  of  this  figure  shows  that  a  strong 
correlation  exists  between  oxidant  concentration  and  temperature.   It  appears  that  patterns 
of  temperature,  which  presumably  reflect  fluctuations  in  radiation  load,  explain  the  oxidant 
variation.   In  addition,  inverse  wind  speed,  since  1965,  also  has  correlated  with  the 
oxidant  fluctuations.   On  the  basis  of  this  figure,  it  would  be  very  difficult  to  demonstrate 
that  the  control  technology  was  dominating  the  observed  oxidant  concentrations.   Figure  2-11 
shows  the  October  analysis  for  carbon  monoxide,  particulate  and  inverse  wind  speed.   Again 
the  correlation  among  these  variables  is  strong  and  no  trend  is  apparent. 

This  elementary  statistical  analysis  illustrates  the  fundamental  difficulty  of  demon- 
strating improvement  in  air  quality  due  to  emission  control.   As  can  be  seen  from  the 
October  analysis,  the  most  important  meteorological  variables,  such  as  inverse  wind  speed 
and  temperature,  dominate  the  pollution  data  and  are  highly  variable  from  year  to  year. 
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Consequently,  significant  improvement  in  quantity  of  pollutants  emitted  from  motor  vehicles 
could  easily  be  masked  by  meteorological  variability.   We  can  state  categorically  that  no 
temporal  trend,  either  up  or  down,  can  be  demonstrated  for  the  majority  of  the  pollutants 
found  in  the  air  above  Sacramento.   There  is  indication  that  hydrocarbon  and  NO-  concen- 
trations have  been  decreasing  during  our  study  period.   The  hopeful  sign  is  countered  by  the 
fact  that  the  associated  pollutant,  oxidant,  does  not  show  an  unambiguous  decrease  for  our 
study  period. 
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CHAPTER  3 
DRIFT,  DIFFUSION  AND  FALLOUT  OF  AN  INSTANTANEOUS  LINE  SOURCE 
OF  HEAVY  PARTICLES  IN  THE  ATMOSPHERIC  BOUNDARY  LAYER 

INTRODUCTION 

In  this  chapter  we  shall  discuss  the  development  of  a  numerical  model  of  the  processes 
by  which  a  line  source  of  gaseous  material  or  heavy  particles,  released  in  the  atmospheric 
boundary  layer,  is  diffused  downwind  and,  in  the  case  of  particles,  deposited  at  the  sur- 
face of  the  earth.   We  have  in  mind  two  applications  for  this  model:   (1)  a  cloud  of  pesti- 
cide or  herbicide  released  by  a  low  flying  aircraft  for  agricultural  purposes  and  (2)  gaseous 
pollutants  and  particulate  matter  emitted  by  motor  vehicles  moving  on  a  highway.   Both  these 
cases  are  serious  environmental  problems  involving  the  health  and  safety  of  individuals  and 
ecosystems  in  the  area  downwind  of  the  source  of  pollution.   In  the  discussion  which  follows, 
we  shall  emphasize  particulate  pollution. 

This  problem  is  one  of  the  most  difficult  conceivable  in  the  field  of  turbulence  and 
diffusion.   Non-isotropic,  boundary  layer  turbulence  is  known  to  differ  in  important  ways 
(Townsend,  1956)  from  the  case  of  homogenious  turbulence,  for  which  a  satisfactory  theory 
is  available  (Batchelor,  1953).   In  addition,  for  the  two  cases  we  are  considering,  the 
effects  of  unstable  and  stable  stratification  are  known  to  be  important  (Slade,  1968; 
Yates,  et   al,  1964).   When  you  add  to  this  system  the  interaction  between  the  turbulent 
wakes,  formed  behind  the  emitting  bodies,  and  the  boundary  layer  flow  the  problem  becomes 
dishearteningly  complex  and  well  beyond  the  current  state  of  knowledge  in  this  field. 
Consequently,  we  adopted  the  idealization  of  the  line  source  as  one  acceptable  simplifica- 
tion.  It  should  be  noted  that  use  of  this  initial  condition  restricts  the  applicability 
of  the  model  to  distances  well  downwind  from  the  release  point. 

The  development  of  this  model  has  been  guided  primarily  by  data  relating  to  the 
deposition  of  aircraft  released  agricultural  chemicals  very  kindly  supplied  by  Professors 
Akesson  and  Yates  of  the  Department  of  Agricultural  Engineering  at  the  University  of 
California,  Davis.   We  found  the  experimental  information  available  from  the  agricultural 
case  far  exceeded  any  that  we  were  able  to  find  for  the  highway  particulate  problem.   We 
do  not  believe  the  problem  of  the  drift,  diffusion  and  fallout  of  agricultural  chemical 
sprays  differs  in  any  important  way  from  the  same  processes  affecting  particulate  matter 


103 


emitted  from  highways.   If  anything,  the  agricultural  chemical  problem  is  more  complex  due 
to  the  possible  importance  of  evaporation  in  changing  the  size  of  the  particles  with  time. 

In  this  chapter  we  shall  develop  the  model  equations  from  the  diffusion  equations  and 
micrometeorological  theory.   Our  basic  objective  has  been  a  model  which  will  predict  depo- 
sition of  particulate  information  from  detailed  knowledge  of  the  emitted  particulate  amount 
and  size  distribution  and  micrometeorological  information.   We  feel  that  the  successful 
formulation  and  testing  of  such  a  research  model  is  the  first  step  towards  developing  an 
operational  model  which  would  be  capable  of  making  useful  estimates  of  deposition  from 
commonly  available  source  strength  and  meteorological  information. 

We  have  deliberately  avoided  the  use  of  arbitrary  and  adjustable  parameters.   We  prefer 
to  use  a  two  stage  development  in  which  we  first  formulate  and  test  a  model  entirely  from 
known  micrometeorological  principles.   We  will  not  introduce  any  corrections  until  the 
basic  behavior  of  the  model  is  revealed.   It  should  be  anticipated,  however,  that  for 
operational  use,  corrections  of  some  kind  would  be  required  by  any  model  due  to  the  great 
complexity  of  the  system  under  consideration. 
Discussion 

The  deposition  of  a  cloud  of  particles  released  in  the  atmospheric  boundary  layer 
can  be  thought  of  as  the  net  result  of  three  processes.   These  are:   gravitational  settling, 
drift,  and  diffusion.   The  first,  gravitational  settling,  is  best  known  and  easiest  to 
calculate.   In  the  size  range  of  interest  here,  diameters  from  1  to  100  microns  (calcula- 
tions for  larger  diameters  can,  of  course,  be  made,  but  these  larger  drops  settle  so  rapidly 
that  the  line  source  approximation  becomes  very  dubious),  the  Stokes  formula  is  valid  only 
for  the  smallest  particles  (Goldstein,  1938).  However,  empirical  information  is  plenti- 
fully available  for  the  entire  range  (Perry,  1950)  so  we  may  regard  the  settling  velocity 
itself  as  well-known  given  reliable  data  on  the  size  distribution  of  the  particles.   Such 
information  is  only  now  becoming  available  for  the  highway  case  (Cahill,  1973).   In  con- 
trast, detailed  size  distributions  have  been  measured  for  some  years  by  Akesson  and  Yates' 
group.   As  we  shall  see,  size  information  is  essential  to  calculations  such  as  we  attempt 
here. 

It  should  also  be  noted  that  the  possible  effect  of  evaporation  on  droplet  size  for 
the  agricultural  case  remains  a  problem  that  we  do  not  consider  in  this  chapter. 
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Even  without  the  formidable  effect  of  turbulent  diffusion,  very  considerable  dispersal 
in  deposition  can  be  achieved  purely  through  the  effect  of  the  differential  fallout  associ- 
ated with  a  size  spectrum.   This  effect  is  illustrated  below  in  Table  1  where  we  list  the 
terminal  velocities  for  particles  of  sizes  ranging  from  10  to  800u  and  deposition  point 
downwind  where  such  a  drop  would  land  if  it  were  released  at  a  height  of  2  meters  with  an 
average  wind,  from  the  surface  to  two  meters,  2  ms   .   As  can  be  seen  the  dispersal  of 
the  smaller  drops  is  quite  large  with  the.  lOu  particles  traveling  over  a  kilometer  before 
depositing. 

Table  1 


Particle  Diameter 

Terminal  Veloc: 

Lty 

Deposition  Point 

(u) 

(cm  s   ) 

(m) 

10 

0.3 

1333 

20 

1.3 

308 

100 

32.0 

12.5 

200 

136.0 

2.94 

800 

340 

1.18 

It  is  also  apparent  why  we  do  not  wish  to  consider  particles  above  100  .   These  sizes  fall 
out  so  fast  that  they  remain  in  the  zone  of  influence  of  the  turbulent  wakes  behind  the 
emitting  vehicles. 

The  second  process,  drift  with  the  mean  wind  depends  critically  on  the  shear  of  the 
wind  with  height  near  the  surface.   This  is  a  phenomena  extensively  studied  in  micrometeo- 
rology  (Goddard,  1971;  Morgan  e_t  al_,  1971)  and  we  can  choose  from  several  alternative 
formulations.   Most  generally, 


az  -  kz  H"/L) 


(l) 


where  U  is  the  wind  speed,  Z  is  the  vertical  ordinate,  U*  is  called  the  friction  velocity, 

t    1/2 

U*   =  (  o/p)    with  t  the  stress  exerted  by  the  wind  on  the  surface,  p  is  air  density, 

k  is  the  von  Karman  constant  with  a  value  of  0.4,  L  is  stability  parameter, 


L  =  - 


PC  TU* 

P 

kgll 


(2) 
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with  C   the  specific  heat  capacity  at  constant  pressure,  T  air  temperature  in  °K,  g 
the  acceleration  of  gravity  and  H  the  flux  of  sensible  heat  transferred  upwards  from  the 
surface  by  turbulence,  and  the  overbar  indicates  a  time  average  of  the  order  of  a  half- 
hour.   The  function  (j>(Z/L)  is  determined  empirically,  and  various  versions  are  available. 
We  have  chosen  to  use  the  "log  +  linear"  approximation  as  the  simplest  which  gives  a 
realistic  effect  of  stability  on  shear.  That  is, 


«j,(Z/L)  =  1  +  a  |   ,  (3) 


where  a  is  a  constant  whose  value  has  been  discussed  by  Priestley  (1963) 
We  are  aware  that  other  approximations  to  the  "phi"  functions  are  probably  more  accurate 
over  a  greater  range  of  stabilities.  At  the  current  stage  in  the  development  of  the 
model,  however,  it  does  not  seem  appropriate  to  introduce  additional  complexity  in  any 
individual  component  until  the  properties  of  the  entire  system  are  understood.   In  any 
event,  other  versions  of  the  "phi"  function  can  easily  be  introduced  when  it  becomes 
apparent  that  additional  accuracy  in  specifying  the  shear  is  warranted. 

The  mean  wind  speed  between  a  release  height,  H,  and  a  reference  height  Z  ,  can 
easily  be  calculated  from  equations  (1)  and  (3).   Z  has  the  significance  of  being  the 
height  where  the  mean  wind  goes  to  zero  according  to  equation  (1).   Empirically,  Z 
is  found  to  be  a  strong  function  of  the  height  of  the  predominant  roughness  elements  of 
the  surface  (Sellers,  1965).   Thus,  this  parameter  offers  the  means  of  introducing  an 
important  aspect  of  the  underlying  surface,  its  roughness.   Integration  of  equation  (1) 
yields 


U=^ 


o 


(4) 


where  U*  and  L  are  treated  as  constants.   It  has  been  found  that  for  large  canopies, 
equation  (A)  describes  the  mean  wind  profile  accurately  only  if  a  "zero  plane  displace- 
ment", d,  is  used  to  shift  the  origin  upwards  to  a  height  comparable  to  that  of  the  canopy 
thickness.   The  displacement  height,  like  the  roughness  parameter,  Z  ,  is  empirically 
determined.   It  should  be  understood,  therefore,  in  equation  (4)  that  the  origin  of  the 
vertical  coordinate  is  at  d.   We  define  the  mean  wind  speed  between  H  as  the  catch  height, 
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Z  ,  as 
c 


U  = 


H  -  Z 


udZ 


(5) 


Evaluation  of  this  integral  with  equation  (4)  substituted  for  the  wind  speed  yields 


U  =  U. 


(H  £nH  -  Z  £n  Z  ) 

7- — ^-—^ -  +^-(H  +  Z   -  2Z  )  -  £n  Z  -1 

(H  -  Z  )       2L       c     o'  o 

c     

£n  -*  +  f  (Z  -  Z  ) 

Z     L    w    o 
o 


(6) 


where  U  is  the  wind  speed  measured  at  the  height  Z  .   Thus,  from  equation  (6)  the  vertically 
w  w 

averaged  wind  speed  can  be  calculated  given  a  wind  speed  measured  at  a  single  height,  the 

surface  roughness  parameters  Z  and  d,  and  the  stability  parameter,  L. 

o 

The  relative  roles  of  stability  and  surface  roughness  in  determining  the  transport 

speed,  U,  may  be  illustrated  by  a  few  order-of-magnitude  calculations  with  equation  (6). 

However,  some  care  is  required  in  formulating  the  calculation  in  order  to  avoid  misleading 

results.   For  instance,  if  we  specify  that  the  release  height,  H,  is  6  meters,  L  =  °° 

(neutral  conditions),  U  =  2  ms   at  Z  =  2  m,  and  a  =  2,  the  numbers  in  Table  2  are  ob- 

w  w 


tained. 


Table  2 


Z  (cm) 
o 


U  (ms-1) 


(u)"1 


0.1 
1 
10 

ion 


2.026 
2.041 
2.112 
3.318 


0.493 
0.490 
0.473 
0.301 


We  show  the  inverse  of  the  transport  speed  because  pollutant  concentrations  should  be 
proportional  to  this  parameter,  other  factors  being  the  same.   Inspection  of  these  numbers 
would  lead  one  to  the  anomalous  conclusion  that  the  transport  speed  increased  with  rough- 
ness.  This  is  at  variance  with  a  physical  concept  of  surface  roughness  being  a  major 
factor  in  bringing  about  the  retardation  of  flow  in  the  atmospheric  boundary  layer.   The 

reason  for  this  peculiar  result  lies  in  specifying  the  measured  wind,  U  ,  within  the 

w 
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layer  of  interest,  H  -  Z  .   Other  factors  being  equal,  the  speed  would  not  be  the  same  at 

o 

this  height  as  the  boundary  layer  flows  over  surfaces  of  varying  roughness.   To  do  a 

complete  examination  of  this  problem  would  require  a  substantial  detour  into  the  theory 

of  the  planetary  boundary  layer  since  it  is  the  gradient  wind  speed  that  should  be  held 

constant  in  the  calculation.   Since  that  problem  is  well  beyond  the  scope  of  this  chapter 

(discussed  in  Volume  II  of  this  report) ,  we  shall  do  a  simpler  calculation  that  reveals 

better  the  physics  of  the  situation.   The  numbers  shown  in  Table  3  result  from  specifying 

the  release  height,  H,  to  be  2  meters,  while  the  measured  wind  speed,  U  ,  is  held  constant 

w 

at  2  ras   at  a  height  of  6  meters.   This  gives  us  the  problem  of  the  response  of  the  lower 
portion  of  the  boundary  layer  to  a  change  of  roughness  while  the  wind  speed  is  held  con- 
stant above.   Here  we  also  show  an  order-of-magnitude  calculation  as  to  the  effect  of 
stability  on  the  transport  speed  for  various  stabilities.   The  case  L  =  +  10  meters 
corresponds  to  moderately  stable  conditions  while  L  =  -  10  meters  to  moderately  unstable 
conditions.   These  results  are  physically  more  satisfying.   The  transport  speed  decreases 
strongly  with  increasing  roughness.   In  comparison,  it  can  be  seen  that  the  effect  of 
stability  on  the  transport  speed  is  quite  small  particularly  for  values  of  Z  10  cm  and 


smaller. 


Table  3 


z 

o 

L  = 

+ 

10  m  ' 

(cm) 

U  (msT1 

) 

u-1 

0.1 

1.512 

0.661 

1 

1.351 

0.740 

10 

1.068 

0.936 

100 

0.476 

2.099 

L  = 


L  =  -  10 


U 


u"1 


u 


u"1 


1.518  0.658 

1.352  0.739 

1.052  0.951 

0.431  2.319 


1.525  0.655 

1.353  0.739 

1.033  0.967 

0.371  2.693 


The  third  process  to  be  treated  by  the  model,  diffusion,  is  the  most  difficult  of 
all.   It  is  sometimes  treated  in  air  pollution  studies  by  introducing  turbulent  diffusivi- 
ties  formally  analogous  to  molecular  dif fusivities  but  orders-of-magnitude  larger.   We 
shall  also  adopt  this  procedure  but  first  we  shall  discuss  some  of  the  limitations  of  the 
approach. 

One  of  the  fundamental  properties  of  turbulent  diffusion  is  that  the  diffusivity  is  a 
strong  function  of  the  scale  of  the  diffusing  "patch"  of  fluid.   This  is  a  consequence  of 
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the  fact  that  for  geophysical  turbulence  the  maximum  of  the  spectrum  of  turbulence  is  almost 
always  found  at  quite  long  wavelengths.   Consequently,  as  a  patch  of  marked  fluid  grows  by 
diffusion,  it  becomes  subject  to  diffusion  by  more  and  more  energetic  scales  of  motion. 
As  L.  F.  Richardson  anticipated,  the  diffusivity  increases  roughly  as  the  4/3  power  of 
the  dimension  of  patch  of  marked  fluid  in  question.   Thus,  diffusivity  is  not  a  constant 
as  is  most  convenient  mathematically.   For  the  problem  under  consideration,  this  means 
that  the  diffusivity  should  increase  with  time  as  the  cloud  of  particles  moves  downwind 
and  increases  in  size  due  to  turbulent  diffusion.   Under  these  conditions,  the  choice  of 
a  constant  diffusivity  leads  to  fundamental  difficulties. 

Another,  more  operational  criticism  that  can  be  made  of  the  diffusivity  concept  is 
that  ordinarily  no  effort  is  made  to  relate  dif fusivities  used  in  air  pollution  studies 
with  micrometeorological  processes.   For  instance,  surface  roughness  and  stability  are 
known  to  exert  great  influence  on  turbulence  in  the  atmospheric  boundary  layer  yet  these 
quantities  are  only  rarely  taken  into  account  in  choosing  dif fusivities.   One  of  our  ob- 
jectives in  the  model  under  discussion  was  explicitly  to  include  micrometeorological 
parameters,  wherever  possible,  so  that  the  influence  of  variable  terrain,  crops  and 
stability  could  be  taken  into  account  in  an  objective  fashion. 
The  model 

We  shall  consider  the  intial  cloud  to  be  an  instantaneous  line  source  at  height  H. 

The  wind  is  taken  to  be  blowing  perpendicular  to  the  line  of  release  and  is  measured  to  be 

U  at  height  Z  .   The  downwind  coordinate  direction  is  called  X.   The  coordinates  of  the 
w  w 

center  Y,  the  plume  at  time  t,  are  X  and  Z  .   The  plume  is  considered  to  be  infinite  in 

P      P 

the  Y-direction. 

We  assume  that  the  center  of  the  cloud  moves  with  the  speed  of  the  wind  averaged  from 
H  to  the  surface,  U.   Then 


and 


Xp  =  Ut      ,  (7) 


Z  -  H  -  ct     ,  (8) 

P 
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where  c  is  the  terminal  velocity  for  a  specified  particle  size  and  X  =  0  and  Z  =  li  at 

P  P 

t  =  0.  In  order  to  calculate  the  downwind  drift  of  the  remainder  of  the  cloud  after  the 
center  reaches  the  surface  (or  catch  height)  we  shall  continue  to  displace  the  center  at 
speed  U,  "underground". 

Since  the  terminal  velocity,  c,  depends  on  particle  size,  the  manner  the  size 
distribution  is  specified  is  critically  important.   We  have  developed  two  schemes.   In 
one,  we  assume  the  particle  size  spectrum  is  Gaussian  so  that  we  need  only  specify  the 
mean  diameter  and  the  standard  deviation  to  determine  the  entire  distribution.   It  is 
necessary  to  truncate  the  spectrum  at  zero  diameter  to  avoid  "negative"  diameters.   The 
continuous  Gaussian  distribution  is  approximated  by  a  series  of  step  functions,  as  the 
number  of  discrete  steps  is  chosen  according  to  the  degree  of  approximation  desired. 
In  the  second  scheme  for  specifying  the  particle  size  spectrum,  the  distribution  is 
completely  arbitrary.   Any  number  of  discrete  steps  can  be  accommodated.   In  this  way,  If 
measurements  of  the  size  distribution  are  available,  the  observed  spectrum  may  be  approxi- 
mated to  whatever  degree  desired.   Thus,  in  this  model  we  replace  the  "real"  cloud  with  a 
continuous  size  distribution  with  a  series  of  homogenious  clouds  each  falling  with  a 
different  terminal  velocity.   As  we  will  describe  in  the  following  pages,  we  repeat  all 
our  drift,  diffusion  and  fallout  calculations  for  each  cloud  and  then  sum  them  all  to 
obtain  the  total  effect. 

We  assume  that  the  mass  flux  in  the  downwind  direction,  F  ,  can  be  represented  by 


F  =  U  q  -  K  -^  (9) 

x         x  3x 


—  _3 

Where  q  is  the  mass  concentration  in  units  of  g  cm   and  K  is  the  turbulent  diffusivity 

2  -1 
in  the  x-direction  in  units  of  cm  s   .   These  two  terms  represent  the  contribution  to 

the  downwind  flux  of  particles  of  mean  drift  and  diffusion.   The  vertical  mass  flux  is 

given  by 


Fz  =  "  C  <  "  Kz  U    •  (10) 


where  c  is  the  terminal  velocity  of  a  given  size  particle  and  is  taken  as  a  positive 
quantity.   These  two  terms  then  represent  gravitational  settling  and  vertical  diffusion. 
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These  equations  are  formulated  so  as  to  give  a  positive  flux  in  the  x  and  z  directions 
and,  therefore,  a  deposition  is  a  negative  flux.   The  overbar  indicates  that  the  quantity 
q  should  be  thought  of  as  an  averaged  quantity  either  through  executing  a  large  number 
of  trials  or  by  averaging  in  the  y  (crosswind)  direction. 

The  equation  governing  mass  concentration  at  some  downwind  point  is 


3t       3x     3z      .2.2 

3x    3z 


(11) 


where  we  have  assumed 


K  =  K  =  K  =  constant. 
x    y 


(12) 


Solutions  may  be  obtained  for  equation  (11)  using  the  boundary  conditions, 

q  ■*■   0  as  t  -*■  0  for  x,  z  >  0 
q  ■+  0  as  t  -*■  «» 

and  defining  a  "source  strength",  Q,  such  that 


Q  = 


-K*> 


q  dxdz     , 


(13) 


i.e.  Q  is  the  total  amount  of  material  released  in  units  of  grams  per  unit  length  along 
the  release  line. 

By  direct  substitution  into  equation  (11)  it  may  be  verified  that  the  solution  is 

I  l 

-|(x  -  Ut)2  +  (z  -  II  +  ct)2] 


-  _2_ 


q  = 


4TTKt 


exp 


4Kt 


(14) 


Differentiation  of  q  with  respect  to  z  and  x  and  use  of  equations  (9)  and  (10)  gives 


F      =   0 
z 


z   -   H  -   ct 
8irKt2 


exp 


-    j(x  -  Ut)2  +   (z  -  H  +  ct)2| 


4Kt 


(15) 
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Fx  =  Q 


x  +  Ut 

8TTKt 


exp 


(x  -  Ut)2  +  (z  -  11  +  ct)2] 


4Kt 


(16) 


The  total  mass  deposition,  per  unit  area  and  normalized  with  the  source  strength,  at  catch 
height,  Z  ,  and  distance  downwind,  x,  is  given  by 


f  ) 

Z   -  H  -  ct\       I  -  [(x  -  Ut)2  +  (Z  -  H  +  ct)2]  I 

c8,Kt>     )  exp|     sr — —    dt 


(17) 


which  we  will  call  the  normalized  deposition. 

It  should  be  noted  that  we  are  identifying  the  total  downward  flux,  calculated  with 
equation  (15),  with  deposition  at  the  catch  height.   This  is  in  marked  contrast  with  the 
usual  procedure  in  which  the  air  concentration  is  calculated  from  equation  (14),  or  its 
equivalent,  at  the  catch  height.   This  number  is  multiplied  with  a  "deposition  velocity", 
an  additional  empirical  parameter  which  supposedly  takes  into  account  the  properties  of 
various  surfaces  in  regard  to  impact  geometry  and  other  processes  of  deposition.   Our 
general  strategy  has  been  to  avoid  the  introduction  of  free  parameters.   Therefore,  we  have 
chosen  not  to  use  the  deposition  velocity  concept. 

The  total  transport,  or  drift,  past  point  x  is  given  by 


x  +  Ut 
8TrKt2 


exp 


I    ~     l(x  -  Ut)2  +   (z  -  H  +  ct)2]  I 
4Kt 


J 


dtdz 


(18) 


where  the  normalized  drift  is  given  in  units  of  mass  per  unit  length  in  the  y-direction. 

We  have  developed  a  numerical  technique  for  integrating  these  integrals.   However, 
the  algorithm  was  so  consuming  of  machine  time  that  we  found  the  economics  of  systemati- 
cally investigating  the  properties  of  the  model  prohibitive.   Consequently,  we  developed 
the  following  analytic  scheme  for  integrating  equation  (17).   The  right  side  of  equation 
(17)  is  of  the  form 


^dt  = 


2    2 
B  +  C 


1  ,   2     -4Kt 

T +  ~  e 
t  / 


dt 


(19) 


where  C  =  C/8ttK   C  -  H  -  Zc/8iTk   B  =  x  -  Ut  and  C  =  Z  -  H  +  ct.   The  integral  may 
be  further  transformed, 
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F 

Q"dt  = 


/c   c  \ 

_i  +  _2 
,  t     2 

\  t  / 


c3  c5t 


dt 


(20) 


where  Cn  =  (Z2  +  H2  +  X2  -  2  Z  H) /4K   Cr  =  (U2  +  c2)/4K,  C£  =  e  4  and 
3     c  c         5  6 

C,  =  [Ux  -  c(Z  -  H)]/2K.   From  the  Bessel  function  literature  (Dwight,  1955),  we  have 
the  relation 


»/. 


:V_1  exp  (-0  XP  -  pX_P)  dx  =  -  (|)   2P  Bv/   (2/F  ) 


(21) 


for  (J  <  0  and  y  <  0,  where  B  (Y)  is  a  modified  Bessel  function  of  Y  of  order  n.   Our  in- 
tegral may  be  written  as 


dt  = 


+  C6C2 


ci    S    -r-c^ 

C6  ("T  +  -f    e  dt  =  C6C1 


t"1"1  exp(-C  t1  -  C.t"1)  dt   , 


0-1    ,   „      1      -I. 
t    exp(-C  t  -  C  t   )  dt 


(22) 


which  are  two  integrals  of  the  proper  form.   We  can  now  write 


_1 
^dt  =  2C6ClBo    (2^S>   +  2C6C2  Vt\     2     \   (2y/SS}        • 


(23) 


where 


2  C6C1  =  2   C1C6  =   4^  eXp 


c(Z      -  H)    -  xU 
c 


-2K 


2   C6C1  "   2   Cl£     =  ^^     6XP 


H   -   Z 
4ttK 


c(Z      -  11)    -  xU 
c 

-2K 


2  ^fs-k 


2  2  2  2—2 

(Z     +  H     +  X     -  2  HZ)    (c     +  U   ) 
c 


and 


C7C:     =    (Z2  +  H2  +  X2  -  2   HZ)/(c2  +  U2) 
3     5  c 
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The  final  form,  then,  which  we  now  use  in  the  model  is 

r 


—  dt  =  (4ttK)     exp 


C(Z   -  H)  -  xU 
c 

-2K 


C  B 


1   7  2    2    2  2—2 

-£     /(Z  +  H  +  X  -  2HZ  )(c  +  U 
2K.     c  c 


+  (H  -  Z  ) 
c 


2   2    2 
ZtH  +  X  -  2HZ 
c  c 


2   —2 

c  +  U 


+  B, 


-r£  /(Z2  +  H2  +  X2  -  2HZ  )(c2  +  U2) 
2k     c  c 


(24) 


We  have  incorporated  tables  from  Dwight  (1955)  for  the  modified  Bessel  functions  of 
order  0  and  1  into  our  algorithm  for  calculating  normalized  deposition.   Excellent 
agreement  was  achieved  between  the  numerical  and  analytic  integration  schemes.   The 
analytic  algorithm  was  found  to  be  two  orders  of  magnitude  faster  than  the  numerical 
scheme. 

The  last  important  parameter  to  specify  in  the  model  is  the  diffusivity  which  we  have 
assumed  constant  in  order  to  be  able  to  use  an  analytic  solution.   We  have  chosen  to  assume 
that  the  diffusivity  for  the  diffusivity  for  the  small  particle  we  are  considering  is  the 
same  as  that  for  momentum.   This  cannot  strictly  be  true  but  we  have  adopted  the  assumption 
as  a  useful  first  approximation.   From  equation  (1),  we  have 


K  = 


kU*Z 

♦(z/L) 


(25) 


We  choose  to  evaluate  the  diffusivity  at  (H  +  Z  )/2.   Using  equation  (4)  to  specify  the 
friction  velocity,  we  obtain 


K.  - 


k  U   (H  +  Z  ) 
w       c 


In  -±     +  f   (Z   -  Z  ) 
Z     L    w    o 
o 


1  +  2-  (h  +  Z  ) 


(26) 


In  summary,  the  following  information  must  be  specified  in  order  to  calculate  down- 
wind deposition:   (1)  the  wind  speed,  U  ,  at  a  particular  height,  Z  ,  (2)  the  atmospheric 

w  w 
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stability  either  by  means  of  the  parameter  L  or  some  other  measure  of  stability  which  can 
be  converted  to  L  such  as  the  Richardson  number  or  the  Stability  Ratio,  (3)  the  canopy 
roughness  characteristics  as  specified  by  Z  and  the  zero-plane  displacement  height,  (4) 
the  properties  of  the  cloud  of  particles,  i.e.,  particle  density,  size  distribution,  and 
release  height,  11,  and  (5)  the  catch  height,  Z  .   It  should  be  evident  that,  in  so  far  as 
the  state  of  the  art  allows,  we  have  successfully  avoided  use  of  arbitrarily  determined 
parameters  in  this  model.   This  may  be  regarded  as  an  excessively  purest  strategy  and  our 
results  may  indeed  not  be  as  good  as  those  achieved  with  models  in  which  various  parameters 
are  available,  such  as  diffusivity  or  deposition  velocity,  to  adjust  the  calculation  to 
correspond  with  observation.   We  feel  that  more  may  be  learned  from  discrepancies  than 
from  agreement  in  this  case. 
Analysis  and  preliminary  validation 

In  this  section  we  shall  illustrate  the  use  of  the  model  and  the  kind  of  analysis  which 
may  be  performed  with  the  output  information.   We  do  not  pretend  that  the  analysis  is  com- 
plete.  It  is  our  opinion  that  far  too  little  attention  is  usually  paid  the  analysis  of  the 
properties  of  models,  considered  apart  from  the  accuracy  of  their  predictions.   In  many 
cases,  detailed  studies  such  as  systematic  sensitivity  analyses  yield  information  which 
may  be  more  valuable  than  a  bare  validation  effort,  particularly  considering  the  general 
difficulty  of  securing  definitive  experimental  data. 

At  this  time  we  have  carried  out  one  validation  study.   It  is  helpful  to  discuss  the 
nature  of  the  data  which  Professors  Akesson  and  Yates  have  made  available  to  us.   Their 
experiments  are  aimed  at  understanding  the  nature  of  all  the  processes,  from  the  nozzle 
to  the  plant,  by  which  aircraft-released  agricultural  chemicals  are  deposited  in  a  target 
area  and  dispersed  into  the  environment.   The  objective  is  to  assist  in  the  development 
of  systems  which  minimize  drift  beyond  the  target  area  and  maximize  the  effectiveness  of 
the  application.   In  these  experiments,  the  aircraft  makes  repeated  passes  at  an  altitude 
of  about  10  feet  upwind  of  the  instrumented  area.   The  multiple  passes  greatly  increase 
the  reliability  of  the  experimental  statistics.   The  released  material  is  designed  to 
simulate  the  physical  properties  of  various  pesticides  and  herbicides  and  contains  accu- 
rately known  amounts  of  fluorescent  dye.   Downwind  deposit  is  measured  by  determining  the 
amount  of  fluorescent  dye  deposited  on  mylar  plates,  while  air  concentrations  are  measured 
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by  means  of  high  volume  samplers.   Cascade  impactors  are  used  to  determine  the  size  dis- 
tribution of  the  particles.   The  deposition  measurements  are  accompanied  by  micrometeo- 
rological  measurements  of  the  temperature  and  wind  speed  profile  and  turbulence  level. 
The  Stability  Ratio,  SR,  is  used  as  a  measure  of  atmospheric  stability,  where 

t"  -  Y 

2     1         5 
SR  =  -=— — -  x  10J        ,  (27) 

U 

where  the  mean  temperature,  in  °C,  is  measured  at  32  and  8  feet,  the  mean  wind  speed  is 
measured  at  16  feet  in  units  of  cm  s 
Model  verification 

We  had  hoped  to  include  in  this  chapter  a  comprehensive  analysis  of  the  model's 
behavior  relative  to  observed  data.   Unfortunately,  the  experimental  effort  and  analysis 
of  the  data  is  so  difficult  and  time  consuming  that  we  have  been  able  to  complete  only  one 
case.   The  difficulty  lies  in  measuring  the  size  distribution  of  the  emitted  particles. 
Inasmuch  as  we  have  no  idea  how  representative  this  single  case  is,  we  have  chosen  not  to 
present  the  analysis  at  this  time.   We  would  have  preferred  to  delay  publication  of  this 
report  until  comprehensive  verification  information  was  available.   However,  since  the 
delay  in  publishing  this  report  has  already  been  extensive,  we  have  decided  to  present 
the  model  without  verification  at  this  time. 
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